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Preface 


The interest lies in a dynamic formulation capable of simulating the flight of flexible aircraft 
on a computer. Such simulations can speed up the aircraft design process. Moreover they are 
indispensable to autonomous aerial vehicles, which require autopilots. 

Computer simulation of the flight of flexible aircraft is no trivial matter. Indeed, they require 
a new paradigm, discarding confining assumptions and adopting potent methodology. This work 
develops such a paradigm, which amounts to treating the aircraft as a single system. To this 
end, it integrates into a single mathematical formulation the disciplines pertinent to the flight of 
flexible aircraft, namely, analytical dynamics, structural dynamics, aerodynamics and controls. 
The unified formulation is based on fundamental principles and incorporates in a natural manner 
both rigid body motions of the aircraft as a whole and elastic deformations of the flexible compo- 
nents (fuselage, wing and empennage), as well as the aerodynamic, propulsion, gravity and control 
forces. The aircraft motion is described in terms of three translations (forward motion, sideslip 
and plunge) and three rotations (roll, pitch and yaw) of a reference frame attached to the unde- 
formed fuselage, and acting as aircraft body axes, and elastic displacements of each of the flexible 
components relative to corresponding body axes. The equations of motion are expressed in a form 
ideally suited for efficient computer processing. A perturbation approach permits division of the 
problem into a nonlinear flight dynamics problem for maneuvering quasi-rigid aircraft and a linear 
“extended aeroelasticity” problem for the elastic deformations and perturbations in the rigid body 
translations and rotations, where the solution of the first problem enters as an input into the second 
problem. The term “extended aeroservoelasticity” refers to a family of problems, each problem 
characterized by an input from a different aircraft maneuver, with the corresponding equations in- 
volving not only both elastic and rigid body variables but also coefficients depending on any given 
maneuver, and hence coefficients depending generally on time. This is materially different from 
the common aeroservoelasticity, which involves for the most part elastic variables alone and is not 
subject to inputs from aircraft maneuvers. The control forces for the flight dynamics problem are 
obtained by an “inverse” process. On the other hand, the feedback control forces for the extended 
aeroelasticity problem are derived by means of LQG theory. A numerical example presents time 
simulations of rigid body perturbations and elastic deformations about 1) a steady level flight and 
2) a steady level turn maneuver. 

It should be pointed out that sufficient details are provided so as to permit interested parties to 
replicate the results. 
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1. Introduction and Literature Review 

This investigation is concerned with a dynamic formulation capable of simulating the flight of 
flexible aircraft. It integrates seamlessly in a single and consistent mathematical formulation all 
the necessary material from the pertinent disciplines, namely, analytical dynamics, structural dy- 
namics, aerodynamics and controls. The unified formulation is based on fundamental principles 
and incorporates in a natural manner both the rigid body motions and the elastic deformations, 
and the couplings thereof, as well as the aerodynamic, propulsion, control and gravity forces. In 
essence, the formulation not only unifies flight dynamics and aeroelasticity, traditionally treated as 
separate disciplines, but, going beyond that, it permits a computer simulation of the response of 
flying flexible aircraft to external stimuli. 

In describing the motion of rigid bodies in space, it is convenient to attach a set of axes to 
the body, where the axes are commonly known as body axes, and express the motions in terms 
of components along these body axes. It is quite common to describe the motion of rigid bodies 
in terms of the translation of the origin of the body axes and the rotation of the body axes; the 
corresponding variables, particularly the rotations, are referred to as quasi-coordinates. If the origin 
of this reference frame coincides with the mass center of the body, the translations and rotations 
are independent of each other. Moreover, if the body axes themselves coincide with the principal 
axes of the body, then the products of inertia are zero, so that the inertia matrix is diagonal. 

The situation is more complicated for flexible bodies, in which case there are basically two types 
of reference frames: 

i. Fixed in the undeformed body - In this case, it is convenient to define the translation of 
the origin of the reference frame and the rotation of the reference frame as the rigid body 
translation and rotation of the body, and regard any displacement relative to the reference 
frame as elastic deformation. 

ii. Moving relative to the undeformed body - In this case, it is common to choose the reference 
axes so that the linear momentum and angular momentum vectors due to elastic deformations 
vanish; axes satisfying these conditions are called mean axes. Because the elastic defor- 
mations depend in general on time, the mean axes are continuously moving relative to the 
undeformed body. Of particular interest are mean axes with the origin at the system mass 
center, because in this case the three types of motion, namely, the translations of the refer- 
ence frame, the rotations of the reference frame and the deformations measured relative to this 
frame, are all inertially decoupled. The price to be paid for the use of mean axes is very steep, 
however, as the constraints defining these axes are not easy to enforce. For this reason, it is 
common to invoke the use of mean axes to justify inertial decoupling, but without enforcing 
the constraints. 

Whereas proper use of mean axes for flexible bodies in vacuum can produce inertial decoupling, 
in the case of aircraft any such benefits are questionable, because the equations of motion remain 
coupled through the aerodynamic forces. Moreover, if one insists on the use of mean axes, then 
the aerodynamic forces must also be expressed in terms of components along the same mean axes, 
which is a very tedious task at best. 

The motion of force-free rigid bodies in space is unstable. Under certain circumstances, the 
motion of rigid bodies can be stabilized by imparting them some spin about the axis of either 
the minimum or the maximum moment of inertia. On the other hand, flexible bodies cannot be 
stabilized if the spin is about the axis of minimum moment of inertia. The preceding statements 
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imply that the motion takes place in vacuum, such as in the case of a spacecraft. Matters are 
entirely different for flexible aircraft, which are neither force-free, nor do they operate in vacuum. 
In fact, aircraft are subjected to aerodynamic and wind forces, and any stabilization is done by 
active means, namely, by the engine throttle and control surfaces, where the latter consist of the 
aileron, elevator and rudder. 

The flight of aircraft tends to be quite diverse, ranging from steady level cruise to complex 
maneuvers, and stabilization requires the use of controls permitting the aircraft to carry out the 
intended maneuver while suppressing any deviations from it, whether in the form of rigid body 
displacements or elastic deformations. Both flight dynamics and aeroelasticity are concerned with 
aircraft stability. But, whereas flight dynamics is concerned mainly with rigid body motions, aeroe- 
lasticity is concerned with vibration and flutter. 

A stability analysis tends to be limited in scope, in the sense that it can only yield a qualitative 
statement about the nature of motion in the neighborhood of an equilibrium state of a system. 
More specifically, it can state whether the equilibrium is merely stable, asymptotically stable, or 
unstable. The time plays no role in a stability analysis. In fact, stability analyses tend to be limited 
to cases in which the time variable can be eliminated, such as when the equations of motion can 
be reduced to an eigenvalue problem. The stability of time-dependent maneuvers can only be 
evaluated numerically. 

To obtain information going beyond stability statements, such as the time response of aircraft 
to external stimuli, it is necessary to undertake a simulation of the equations of motion, which 
amounts to the integration of the equations of motion. Such a dynamic simulation can be used for 
parametric studies in preliminary design. Moreover, it can be used to determine aircraft perfor- 
mance, thus reducing the time required for flight testing by “flying the aircraft on a computer.” 

The choice in this paper is to work with a reference frame attached to the undeformed aircraft, 
which has many advantages over mean axes. But, because the elastic deformations prevent the 
origin of a frame attached to the undeformed body from coinciding with the mass center and the 
axes themselves from coinciding with the principal axes for all times, there is no preferred choice 
of a reference frame; we base the choice on geometric considerations. In particular, we attach a 
set of body axes to the undeformed fuselage, where one of the axes is along the symmetry axis. 
For convenience, sets of body axes are also attached to the other flexible components, such as the 
wing and the empennage. Ultimately, however, all motions are referred to the fuselage body axes, 
which act as a reference frame for the whole aircraft. 

The mathematical formulation is based on equations of motion in terms of quasi-coordinates 
derived earlier by the first author for flexible spacecraft and later adapted by him to flexible aircraft. 
The formulation is hybrid in nature, in the sense that it consists of ordinary differential equations 
for the rigid body translations and rotations of the aircraft as a whole and boundary-value problems 
for the elastic deformations of the flexible components of the aircraft, namely, the fuselage, wing 
and empennage. For practical reasons, the distributed variables describing the boundary-value 
problems for the individual components are discretized in space, obtaining a relatively large set 
of second-order ordinary differential equations for the whole aircraft. The discretization process 
amounts to representing the distributed variables as finite series of known space-dependent shape 
functions multiplied by time-dependent generalized coordinates. The derivation of the equations 
of motion in conjunction with the extended Hamilton principle requires expressions for the ki- 
netic energy, potential energy and virtual work, all scalar quantities. In turn, the kinetic energy 
requires the velocity of every point of the aircraft, which can be obtained by means of an orderly 
kinematical synthesis, going from the fuselage to the wing and to the empennage. The potential 
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energy is merely equal to the strain energy. Moreover, the aerodynamic, propulsion, control and 
gravity forces are accounted for through the virtual work. Rather than deriving first hybrid equa- 
tions of motion and then discretizing them in space, it is perhaps more expeditious to carry out the 
discretization directly in the kinetic energy, potential energy and virtual work, thus obtaining the 
desired set of ordinary differential equations for the whole flexible aircraft without the need to de- 
rive boundary- value problems. For integration of the differential equations and for control design, 
it is necessary to transform the set of second-order differential equations into a set of first-order 
differential equations, namely, into a set of state equations. It turns out that, for the problem at 
hand, it is more convenient to work with momenta rather than with velocities. Although the result- 
ing first-order equations actually represent phase equations, we shall continue to refer to them as 
state equations. 

The simulation of the flight of an aircraft amounts essentially to integration of the state equa- 
tions. Because of various nonlinearities involved, such as those due to rigid body motions and 
aerodynamic forces, the integration must be carried out numerically on a computer. In one way or 
another, computer integration can only be done in discrete time, which raises the question of the 
size of the sampling period, or time step. Of course, the answer depends on the desired accuracy 
of the simulation, and it is intimately related to the dynamic characteristics of the system. If the 
aircraft is to be controlled by an autopilot, then the simulation must be carried out in real time. If 
the dynamic characteristics are such that the time step must be relatively short, perhaps of the order 
of 0.01s, most aerodynamic theories in current use must be ruled out, as the computation of the 
dynamic pressure over the entire aircraft is sure to take considerably longer time than that. Hence, 
a new method for computing the dynamic pressure must be developed, one characterized by high 
computational speed, even if some accuracy must be sacrificed. Moreover, the method for com- 
puting the dynamic pressure must be compatible with the method for modeling the airframe. If the 
formulation is to be used for aircraft design, then real-time simulation may not really be necessary, 
although on-line simulation may. But, the size of the sampling period, which is determined by the 
system dynamic characteristics, remains the same regardless of whether the simulation is in real 
time or only on-line. The implication is that the computation of the dynamic pressure must still be 
relatively fast. Indeed, a mere 10 s simulation requires 1,000 time steps. Hence if the computation 
of the dynamic pressure takes one hour, the simulation requires over 40 days. This demonstrates 
the need for a method for computing the dynamic pressure in a very short time period. In this 
regard, a reasonably accurate approximate method may be acceptable. 

As indicated above, the equations of motion for a flying flexible aircraft are nonlinear, where the 
nonlinearity is due to the rigid body motions and the aerodynamic forces. Moreover, the equations 
tend to be of high order, the order depending on the discretization procedure employed. Hence, one 
can expect difficulties both with a stability analysis and with control design. In addition, difficulties 
can be experienced in the integration process, because some of the variables describing the aircraft 
rigid body motions tend to be large and the variables describing the elastic displacements tend to 
be small. A perturbation approach to the solution can obviate many of these difficulties. More 
specifically, the solution can be represented as the sum of a zero-order part for the large rigid body 
variables and a first-order part for the small elastic variables and perturbations in the rigid body 
variables, where the zero-order quantities are larger than the first-order quantities by at least one 
order of magnitude. Then, the equations of motion can be separated into a zero-order problem 
for the rigid body motions alone and a first-order problem for the elastic displacements and the 
perturbations in the rigid body motions. The state equations for the zero-order problem are of 
order 12 at most; they can be identified as the equations of flight dynamics and can be used to 
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describe aircraft maneuvers. On the other hand, the state equations for the first-order problem 
are of order 12 + 2 n e , where n e is the number of elastic degrees of freedom; they represent the 
extended aeroelasticity equations, where “extended” is in the sense that they include not only 
the elastic displacements but also perturbations in the rigid body variables, where the latter are 
sometimes referred to as “body freedoms.” 

The flight dynamics equations are in general nonlinear and describe the translations and rota- 
tions of the aircraft as if it were rigid. They can be used to design given maneuvers of an aircraft, 
which amounts to solving an “inverse” problem. In the commonly encountered direct problems 
in dynamics of rigid bodies, the forces are given and the equations of motion are solved for the 
state, i.e., for the positions and velocities. In the context of the present formulation unifying flight 
dynamics and aeroelasticity, however, the state representing a desired maneuver is given and the 
problem amounts to determining the engine thrust and the control surface forces permitting the 
realization of the maneuver; this represents an inverse problem. On the other hand, the extended 
aeroelasticity equations are linear, but they contain the state and forces from the flight dynamics 
problem as coefficients and as an input. Hence, there is a set of extended aeroelasticity equations 
for every conceivable aircraft maneuver. If the flight dynamics problem represents steady level 
cruise or a steady level turn maneuver, then the zero-order state and forces are constant and the 
system of extended aeroelasticity equations is linear time-invariant. In this case, the state equa- 
tions lend themselves to a standard stability check, such as one based on the roots of the eigenvalue 
problem, to control design by commonly used techniques, such as the LQG method, and to ready 
integration for simulation of the aircraft response to external stimuli. If the flight dynamics prob- 
lem represents a time-dependent maneuver, such as the transition from one steady state to another, 
then the zero-order state and forces depend on time and the extended aeroelasticity state equations 
are linear time- varying, which precludes a standard stability analysis. However, the state equations 
still permit control design and response simulation. 

The following literature review should help relate the present paper to previous investigations: 
Although flight dynamics and aeroelasticity have been developed traditionally as separate disci- 
plines, the need for considering interacting efforts was recognized quite early. 1 ' 3 Still, relatively 
few attempts have been made to link the two disciplines, and when such attempts were made almost 
invariably the scope was quite limited. This lack of interest in linking flight dynamics and aeroe- 
lasticity can be attributed to a reluctance to increase the complexity of the problem to a significant 
extent at a time when powerful computers capable of solving such problems were not available. 
As a result, problems combining flight dynamics and aeroelasticity effects have tended to be sub- 
jected to many simplifying assumptions designed to permit largely analytical solutions. In one of 
the first references on the subject, Bisplinghoff and Ashley 4 derived scalar equations of motion for 
an unrestrained flexible vehicle. The equations consisted of three inertially decoupled sets, one for 
the rigid body translations, one for the rigid body rotations and one for the elastic deformations, 
the latter expressed in terms of aircraft structural natural modes. Although not stated explicitly, 
this implies the use of principal mean axes with the origin at the vehicle mass center. Moreover, 
the inertia matrix was assumed to be constant, which implies that the contributions from the elastic 
deformations to the inertia matrix were ignored. Aerodynamic forces for the case of small dis- 
turbed motions from steady rectilinear flight were given in terms of an influence function for an 
unrestrained aircraft. An integrated analytical treatment of the equilibrium and stability of flexi- 
ble aircraft was presented by Milne. 5 In Part I, he derived linearized equations of motion about a 
steady state, assuming not only that the elastic deformations but also the rigid body translations 
and rotations were small. Although the constraint equations defining the mean axes were given, 
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the formulation seems to have used body axes attached to the undeformed aircraft. The equations 
are expressed in a vector-dyadic form that does not permit a ready check for missing terms and, 
more importantly, does not lend itself to ready computer solutions. In Part II, the general analysis 
was applied to the study of equilibrium and longitudinal stability about equilibrium of an aircraft 
having longitudinal flexibility only. A monograph by Taylor and Woodcock 6 consists of two parts 
representing different approaches to the same problem. In Part I, Taylor presents a very lucid sum- 
mary of the equations of motion for deformable aircraft derived by Bisplinghoff and Ashley 4 and 
by Milne. 5 Following a reduction to scalar form, the equations are simplified to permit the study 
of some special cases. In Part II, Woodcock uses an unorthodox form of Lagrange's equations 
to derive scalar perturbation equations of motion about a given “datum motion,” not necessarily 
corresponding to steady level rectilinear flight; the equations are in terms of body-fixed axes. The 
question of aerodynamics receives scant attention in both parts. An extensive report by Dusto et 
al., 7 resulting in a computer program known as FLEXSTAB, integrates flexible body mechanics 
with a low frequency aerodynamics employing linear influence coefficients. The flexible aircraft 
mechanics uses free vibration modes superimposed on rigid body dynamics. Aerodynamic influ- 
ence coefficients are derived using a paneling scheme lending itself to empirical corrections. The 
equations are expressed in terms of steady perturbations about a reference motion to determine dy- 
namic stability by characteristic roots or by time histories following an initial perturbation or some 
gust disturbance. There are two major concerns. The first consists of the fact that the structural 
dynamics formulation is in terms of mean axes and the aerodynamics is in terms of a different 
set of axes, namely, “fluid axes,” where the latter move with a steady velocity relative to the un- 
deformed aircraft body axes; using two different sets of axes in the same formulation, without 
making the necessary transformation from one set to another, is a very questionable proposition. 
The second source of concern is the time required to run FLEXSTAB (see Ref. 24). Several ana- 
lytical methods for mathematical modeling of aircraft active control system design are described 
by Roger, 8 placing the emphasis on aerodynamics. Inconsistencies in control configured vehicles 
are highlighted by Schwanz. 9 He suggests that familiarity of flight control specialists with a broad 
spectrum of pertinent disciplines, including aerodynamics, structures, modem dynamics and con- 
trol, can minimize and perhaps avoid altogether these inconsistencies. Free-free dynamic analyses 
of forward swept wing aircraft by Miller, Wykes and Brosnan 10 have shown that the static aeroe- 
lastic divergence exhibited by a cantilevered forward swept wing is replaced by a low-frequency 
flutter mode due to coupling between the wing divergent mode and the aircraft short-period mode. 
This coupling is shown to have detrimental effects on flying qualities, ride qualities and gust loads, 
but these effects can be minimized by an active flutter control system. In the same spirit, Weisshaar 
and Zeiler 11 discuss the importance of including aircraft rigid-body modes in the aeroelastic anal- 
ysis of forward swept wing aircraft. They show that body-freedom flutter and aircraft aeroelastic 
divergence, not wing divergence, are the primary aeroelastic instabilities. Rodden and Love 12 point 
out that equations of motion derived using mean axes for the inertial terms and axes attached to the 
undeformed structure for the flexibility terms are likely to be incorrect; such flexibility terms are 
obtained when using structural influence coefficients. Cerra, Calico and Noll 13 developed a linear 
model of an elastic aircraft providing the capability of analyzing the coupling between the rigid 
body motions and the elastic motions. The model can be used for stability and control analyses. 
As in Ref. 4, the rigid body translations, rigid body rotations and elastic deformations are assumed 
to be inertially decoupled. A framework for constructing simulation models for flexible aircraft is 
described by Arbuckle, Buttrill and Zeiler. 14 The objectives are to increase simulation model fi- 
delity and to reduce the time required for developing and modifying these models. The framework 
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has been applied to the development of an open-loop F/A-18 simulation model. Buttrill, Zeiler 
and Arbuckle 15 considered a mathematical model integrating nonlinear rigid body mechanics and 
linear aeroelasticity in conjunction with Lagrangian mechanics to derive the equations of motion 
for flexible aircraft. Undamped vibration modes satisfying first-order mean axes constraints were 
used as generalized coordinates. Considering a model of an F/A-18 aircraft, elastic modes sig- 
nificantly affected by inertial coupling were found to be aerodynamically decoupled from the rest 
of the model. Zeiler and Buttrill 16 used the extended Hamilton principle to derive equations of 
motion for a flexible body. The equations include inertial terms due to flexibility, as well as terms 
coupling rigid body and flexible momenta. In addition, a nonlinear strain-displacement relation 
permits centrifugal stiffening to be taken into account. The equations are used to simulate the 
motion of a structure spinning initially about an unstable principal axis in gravity-free vacuum. 
Using Lagrange's equations, Waszak and Schmidt 17 derived the equations of motion for a flexible 
aircraft. The strip theory was used to obtain closed-form integral expressions for the generalized 
aerodynamic forces. Moreover, the use of mean axes permitted inertial decoupling of the rigid 
body translations, rigid body rotations and elastic deformations, the latter being expressed in terms 
of aircraft vibration modes. The modeling method was applied to a generic elastic aircraft, and the 
model was used for a parametric study of the flexibility effects. Nonlinear equations of motion for 
elastic panels in an aircraft executing a pull-up maneuver of given velocity and angular velocity 
were derived by Sipcic and Morino. 18 The effect of the maneuver on the flutter speed and on the 
limit cycle amplitude was discussed for various load conditions. Accurate modeling of aeroelas- 
tic vehicles, with emphasis on the rigid body and elastic degrees of freedom, was discussed by 
Waszak, Buttrill and Schmidt. 19 A comparison of the approach of Ref. 17 on the one hand and that 
of Refs. 15 and 16 on the other hand was presented and various model reduction techniques were 
reviewed. An integrated analytical framework for modeling elastic hypersonic flight vehicles was 
developed by Bilimoria and Schmidt. 20 A Lagrangian approach was used to derive equations of 
motion including rigid body motions and elastic deformations, as well as effects due to fluid flow, 
rotating machinery, wind and a spherical rotating Earth. The elastic deformations are represented 
in terms of modal coordinates relative to mean axes. A paper by Olsen 21 reveals new insights in 
the aeroelasticity and flight mechanics of flexible aircraft by obtaining and solving the equations 
of motion for an accelerating, rotating aircraft. General equations in terms of quasi-coordinates are 
first obtained and then reduced to the case of a “flat” airplane. The influence of gusts on the dy- 
namics of large flexible aircraft is analyzed by Teufel, Hanel and Well, 22 who present an integrated 
flight and aeroelastic control law reducing gust sensitivity. Moreover, the control laws, designed by 
^-synthesis, are such that flight maneuvers do not excite elastic motions. Konig and Schuler 23 de- 
scribe how an integral model for large flexible aircraft can be derived and how an integral control, 
covering flight control, load control and structural mode control, can be designed by multiobjec- 
tive parameter optimization. Samareh and Bhatia 24 presented a unified three-dimensional approach 
that reduces the number of interactions among various disciplines by using a computer-aided de- 
sign model. Results were presented for a blended wing body and a high-speed civil transport. 
Schmidt and Raney 25 considered the effects of flexibility on the flight dynamics of large flexible 
aircraft. In particular, when the frequencies of the lower elastic modes approach those of the rigid 
body modes the handling characteristics can suffer and the flight control system design tends to 
become significantly more complex. Expressing the motion in terms of components along mean 
axes, they add the flexible degrees of freedom to an existing simulation model of the vehicle’s rigid 
body dynamics. The NASA Langley Research Center simulation facility was used to obtain the 
dynamic response of two different aircraft. 
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With some exceptions, the equations of motion in Refs. 3-23,25 were derived either by means of 
Newtonian mechanics or by means of standard Lagrange’s equations. These approaches are more 
suitable when the motions are expressed in terms of inertial axes and/or when the rotations are in 
terms of Euler's angles. Yet, in the case of aircraft it is more convenient to express the motion 
in terms of components along body axes. In this regard, we should point out that this is common 
practice in flight dynamics, in which case the angular veolcities in terms of body axes are the well- 
known roll, pitch and yaw. Of course, equations in terms of inertial axes and/or Eulerian angles can 
always be transformed into equations in terms of body axes through coordinate transformations. It 
is appreciably simpler, however, to derive the equations of motion directly in terms of body axes, 
which can be done through the use of Lagrange's equations in terms of quasi-coordinates. 26 

Motivated by problems in dynamics of spacecraft with flexible appendages, Meirovitch and 
Nelson 27 derived for the first time hybrid (ordinary and partial) differential equations of motion 
coupling rigid body rotations and elastic deformations. The elastic deformations were measured 
relative to a set of body axes attached to the undeformed spacecraft and the rotational motions were 
in terms of quasi-coordinates. The explicit formulation of Ref. 27 was extended by Meirovitch 28 
to a generic whole flexible body by deriving a set of hybrid equations of motion in terms of quasi- 
coordinates, treating for the first time translational velocities as quasi-velocities; the equations were 
then cast in state form. The equations of motion in terms of quasi-coordinates of Ref. 27 were used 
by Platus 29 to derive coupled equations of motion governing the aeroelastic stability of spinning 
flexible missiles. The coupling between the elastic deflections and rigid-body motions was nonlin- 
ear, but the equations were linearized so as to permit a stability analysis. The developments of Ref. 
28 were extended by Meirovitch 30 and Meirovitch and Stemple 31 to flexible multibody systems. 
Then, the approach of Refs. 28, 30 and 3 1 was used by Meirovitch 32 to produce a definitive unified 
theory for the flight dynamics and aeroelasticity of whole aircraft. Generic state equations describ- 
ing the flight of flexible aircraft were first derived in hybrid form and subsequently discretized in 
space. Then, using a perturbation approach, the discrete state equations were separated into a set of 
nonlinear flight dynamics equations for the rigid body variables and a set of linear extended aeroe- 
lasticity equations for the elastic variables and perturbations in the rigid body variables. Nydick 
and Friedmann 33 applied the equations of motion in terms of quasi -coordinates derived in Ref. 28 
to the analysis of a hypersonic vehicle in free flight. To this end, they simplified the equations by 
considering only the pitch and plunge rigid body degrees of freedom and small elastic displace- 
ments. The nonlinear equations were linearized about a trim state obtained by using a rigid body 
trim model and steady hypersonic aerodynamics. Flutter derivatives were calculated by means of 
piston theory. The generic formulation of Ref. 32 was used by Meirovitch and Tuzcu 34 to carry out 
the derivation of explicit equations of motion in terms of quasi-coordinates for a flexible aircraft 
model and to cast the equations in a special state form suitable for simulation on a computer. Due 
to relative ease of integration into the unified formulation and computational speed advantages, the 
aerodynamic forces were derived by means of strip theory. Then, equations for flight dynamics 
and extended aeroelasticity were derived. An approach entirely different from that in Ref. 34 is 
proposed by Fomasier et al. 35 Indeed, Ref. 35 is concerned essentially with the fluid-structure in- 
teraction in a flexible aircraft. To this end, it uses “temporal and spatial algorithms” to make two 
independently developed computer codes, one for the aerodynamics (CFD) and one for structural 
mechanics (CSM), work together. The scope of Ref. 35 is relatively limited, as the aircraft is as- 
sumed to follow a known preset trajectory, so that there are no rigid body degrees of freedom, and 
there are no controls. Worthy of notice is the fact that several 5 s simulations, including some of 
the wing tip displacement, took approximately 35 hrs on a 32-processor Sgi computer. 
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The present paper represents an extension of the developments in Ref. 34. In addition to some 
modeling refinements, the paper contains a numerical example for a model of a flexible aircraft 
containing 76 states, 12 rigid body states and 64 elastic states. Two flight dynamics problems are 
considered, the steady level cruise and a steady level turn maneuver. The corresponding extended 
aeroelasticity problems are derived and used to design feedback controls guaranteeing the vanish- 
ing of the rigid body perturbations and the elastic vibration, and hence the stability of the maneuver 
and the comfort of the flight. The control design consists of a linear quadratic regulator in conjunc- 
tion with a stochastic observer. The integrated process is demonstrated by means of a numerical 
example including a variety of rigid body and elastic displacements time simulations together with 
the corresponding controls time histories, all carried out on a 1 GHz PC using MATHEMATICA. 
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2. Hybrid Equations of Motion in Terms of Quasi-Coordinates 

We regard the aircraft model shown in Fig. 1 as a flexible multibody system subjected to gravity, 
aerodynamic, propulsion and control forces, where the bodies can be broadly identified as the fuse- 
lage, wing and empennage. The motion of the aircraft can be conveniently described by attaching 
a reference frame xjyjZj to the undeformed fuselage (Fig. 1), as well as corresponding reference 
frames x w y w z w and x e y e z e to the wing and empennage, where the various reference frames repre- 
sent respective body axes. Then, the motion can be described by six rigid body degrees of freedom 
of the fuselage body axes, three translations and three rotations, and by the elastic deformation of 
every point of each flexible component relative to the respective body axes. 

From Ref. 32, and making provisions for members in torsion, as well as for structural damping, 
the hybrid equations of motion for the whole flexible aircraft in terms of quasi-coordinates have 






dLj dT ui f 

+ ~Z~. b 4tUi — Uj, 



( 1 ) 


+ 


i = / (fuselage), w (wing), e (empennage) 


where 


L = Lagrangian for the whole aircraft 
V/, = vector of translational, angular quasi-velocities of x/y/Zf 

V/,Uf = skew symmetric matrices derived from V f , u>j 

Cf — matrix of direction cosines between x f y f Zf and XYZ (inertial axes) 

R / = R f(Xf, Yf, Zf ) = position vector of origin Of of XfyfZf relative to XYZ 
Ef = matrix relating Eulerian velocities to angular quasi-velocities 
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Of = symbolic vector of Eulerian angles between x/y/Zf and XYZ 
Uj, v, = elastic displacement, velocity vectors for body i 

oti = elastic angular displacement, velocity vectors for body i 
Li = Lagrangian density for body i exclusive of strain energy 
^uii F vi = Rayleigh’s dissipation function 37 densities for body i 
C ui , = matrices of stiffness differential operators for body i 

F, M = resultant of gravity, aerodynamic, propulsion and control force, moment vectors 
acting on the whole aircraft in terms of fuselage body axes components 

Uj, 'J'i = resultant of gravity, aerodynamic, propulsion and control force, moment density 
vectors for body i 


Assuming that axes xjyfZf are obtained from axes XYZ through the sequence of rotations Y 
about Z to axes x iyiZi, 0 about y\ to x 2 y 2 z 2 and <t> about x 2 to Xfy/Zf, the matrices Cf and Ef 
have the form 36 


Cf = 

E f = 


cipcO sxpcd — s 6 

cipsOsY — sipcY sipsOsY + cipc4> cOscf) 
C1pS0C<fr + s^s 0 SIpsOcY — C0S0 C0C(f> 

1 0 -s9 ' 

0 c (f) c9s<}> 

0 — s<j> c6c(f) 


(2) 


in which s = sin, c = cos. The elastic displacement vectors U* and i/>. are subject to given 
boundary conditions at the interface between bodies. Equations (1) involve the Lagrangian L = 
T — V, in which T is the kinetic energy and V the potential energy, the Rayleigh dissipation 
function densities X ui and which contain the information about the structural damping, and 
the stiffness operators C ui and which are related to the strain energy. Before more explicit 
equations of motion can be derived, it is necessary to produce these quantities. The kinetic energy 
for the whole aircraft can be written as 


T = Tf + T W + T e (3) 

where 

Ti = \ fVfVidmi, i = f,w,e (4) 

are kinetic energies of the individual components, in which V, are velocity vectors of typical points 
in the components and drrii are corresponding mass differential elements. The velocity of a point 
in the fuselage can be written as 

V/(r/, t ) = V f (t) + [f f + u f ( r f , t)] T [u f (t) + <2/(17, *)] + v /( r /> *) 

= V f + (ff + Uf) T Vf + rja f + Vf ( 5 ) 

where r/ is the nominal position of the mass element dm /, f/ and Uf are skew symmetric matrices 26 
corresponding to r/ and u / and v/ and a/ are elastic velocities associated with bending and tor- 
sion, respectively. Then, denoting by C w the matrix of direction cosines between x w y w z w and 
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XfVfZf, the velocity of a point on the wing has the expression 

Vu,(r w ,t) = C w V f (r fw ,t) + r w C w [rif(Tf w ,i) + atf(rf w ,t)] + [ t)] T [u> w (t) 

+ a ti) (r UM f)] +v w (r w ,t) 

— Gyj "V f “1“ [ Gyj (x fw ”t"" U'f'w) ""f"" ) Gw ] ^ f t’yj C yy ( f w fw ) 

+ c w ( V fw + rj w a fw ) + r^ctu, + v u , (6) 


in which Tf w is the radius vector from the origin of Xf tjfZf to the origin of x w y w z w , 

ttfw = [0 — du fz /dxf du fy /dxf}\l^ (7) 

is the angular velocity of the fuselage at r/ w due to bending and af W = [af W 0 0] is the elastic 

velocity of the fuselage at r/ u , due to torsion. The velocity V e (r e , t) of a point on the empennage 
can be obtained from Eq. (6) by simply replacing w by e. The total kinetic energy can be expressed 
in the general form 


T = \ JV'fVfdmf + I JV^V w dm w + § JVjV e dm e (8) 

The potential energy can be expressed in terms of the operators Cui, C^i ( i = f, w, e), but is more 
conveniently expressed as the strain energy. Moreover, the Rayleigh dissipation function densities 
can be expressed in the form 


d^ui — o t'ui E !■ i 


(Pu l&Uj 
dx? dx? ' 


Eipi — o G <J i 


*£ diffj 
dxj dxi 


i = f,w,e 


(9) 


where Cui, c ^ are damping functions and EI t , GJi are flexural and torsional rigidities (i = 

f,w,e). 

Equation (8), the potential energy V, the functions E u i and and the operators C ui and C vi , 
when inserted into Eqs. ( 1 ), permit the derivation of explicit hybrid equations of motion. Because 
for all practical purposes it is not feasible to work with hybrid equations, we do not pursue this 
subject any farther, and approximate instead the partial differential equations by sets of ordinary 
differential equations. 
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3. Spatial Discretization of the Distributed Variables 


For the most part, aircraft are modeled as discrete systems, either by regarding the inertia and 
stiffness properties as lumped from the onset or by describing their elastic motions in terms of 
aircraft structural modes. For undamped structures in vacuum, the use of structural modes yields 
complete decoupling, i.e., it yields a set of independent modal equations of motion. In the case 
of aircraft, however, complete decoupling is not possible, even when aircraft structural modes are 
used, as the aerodynamic forces guarantee that the equations of motion remain coupled. Spatial 
discretization using aircraft structural modes not only does not offer any particular advantage but 
also has the disadvantage that some geometric details of the aircraft are lost in the process. Hence, 
a discretization procedure that does not require the structural modes, which may not be readily 
available and/or may not be compatible with the rest of the formulation, and retains some sense of 
the aircraft geometry seems desirable. Consistent with this, we consider spatial discretization of 
the individual aircraft components separately. To this end, we use either the Galerkin method or 
the finite element method 37 and introduce the expansions 


Ui(Ti,t) = $tti(ri)q Mi (t), = ^(r^q^f), i = f,w,e (10) 

where and are matrices of shape functions and q m and are corersponding vectors of 
generalized coordinates. Note that the choice of shape functions is very important. Indeed, for 
accurate modeling, the shape functions must reflect the mass and stiffness characteristics of the 
components. Some guidelines concerning the choice of shape functions can be found in Ref. 37. 
Moreover, we denote the associated generalized velocities by 


S ui(t) = <Mf), S #i(t) = q i = f,w,e 


(H) 


In anticipation of later needs, we write the velocity vectors for points on the individual components 
in the two discrete forms 


V/(r/,f) = V/ + ( ff + $ u f<i U f) T utf + $ u fS u f + 

= "V f + rjuff 4 - u>f$ u fq u f + $ u fS u f + 

V w (r w ,t) = C w V f -b [C w (ff w “I - $ u , fw q U f ) -b (r w -b ^uu'quiu ) Ciu]^/ 

~b "b C W Q U f W )Suj “b &UW&UW 

+ ipfw "b C w r + T w ^ 1 p w S 

= C W V f + {C w f + r w C w }u) f + CwUfQufwQuf ~b uwQuw 

"b if w "b Cyji&ufw^Suf “b ®ULt ! 

"b ~b C w r + T ipwSxliw 


in which 

A = 


0 0 
0 —d/dxf 
d/dxf 0 


> 


ti fy o; 

o T tit. 


) ^ ufw ^ufi.Tfw') 


( 12 ) 


(13) 


Expression analogous to can be written for <J>^/ and $^f w . Moreover, V e (r e t) can be obtained 
from V w (r w t) by simply replacing w by e. Inserting the first forms of Eqs. (12) into Eq. (8) and 
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( 14 ) 


carrying out the indicated operations, we can write the kinetic energy in the compact form 

T = \ fVjV fdrrif + ± fVlV w dm w + ± fVjV e dm e = |V r MV 

where V = [Vj u> T { sL s T uu , s£ e s^ f sj u , sj e ] r = [Vf ... Vj] T is the discrete system velocity 
vector and M = [M^J is the system mass matrix, a matrix that can be expressed in partitioned 
form with the submatrices 

M n = ml, M n = S T , 

M\z J^u/dmf 4~ J C u ,{t w C w AQ U f w -)- C w $ U f w )dm w 4~ J C e (r e C e ^^ u j e 4- C e $ u j e ) d,fn e 
M u = C'lJ $ uw dm w , M15 = Cj f $ ue dm e 

M\q Jfj^^fdTTif 4“ C u , J ( f w C w 4" c w f fw % fw )dm w 4~ C e J (v e C e ^^f e 

+ C e fJ e ^ti>fe)dm e 

M u = Clfrl t $^ w dm w , M l8 = Cj JrJ% € dm e 
M 21 = M^2, M 2 2 = J 

M 2 z = J(?f 4- 4- J[C w (ff w 4- $ u fw < iuf) T 4- (r w 4- ^UW Quit' ) ] 

^ iTw^ w ^^ufw A v^ufw)d J TYl w “1“ f[^e(ff e A ^li/eQu/) (f e “I - ^ueQue ) T CeY 

x (fJC e A$ ufe + C € $ ufe )dm e 

A^24 J \Cyu{v fw “I” ^llfwQuf^ “1“ (r W “I - Cyj\ Quwd'TH'w 

M25 J^C e (v f e “h ^u/eQii/) + (f e + ^tieQue) Ce] ^ued^R e 

A/26 f (j* f A ^ufQuf)f f ^'tpfdTUf + f[C w (ff w + $ u fwQuf) A (T w “1“ *&uwQ.uw) C u ,] 

-X ( w C44; ^ xp fw A f fw ^xftfxv ) dvn w ~\~ i /[C e (ry e -f- $ u j e cj^y) -|- (r e -h $ iteQue fCeY 

X (/*e A C e ry e $^y e )rf?72 e 

A/27 = ^[Cy; {r fw A ^u/ujQ/u/) A ( A ) ^w\ ^ w *&xj>wdTfl w 

M 2S = J[C € (ff e A $ u fe<luf) T + (fe + $ neQuef C e } T f J % € dm e 

A /31 = Myy M 32 = 

A/33 = f $l f $ufdm f A J (rlC w A$ ufw + C w $ ufw ) T (rlC w A<Z> ufw A C w $> ufw )dm w 
A / (fJCeA$ U fe A C e $ U f € ) T (rJC e A& u f e A C e $ u f e )dm e 
A/34 — f (f^C w A^ u f w + C w Q u f W ) T Q uw dm w ^ 

A/35 = J (rJC e A<P U f e + C e $ u f e ) T Q ue dm e 

A/36 J^ufff %fdmj + f (f w C' W A*P u fw A C w Q u f u ?) (v w C w ^ x pf w -f- Cwf fw*&xpfw )dTn w 

+ /(r^C e A$ w/e + C € ^ uf€ ) T (rJC € ^ f€ + C e r T }e % fe )dm e 
A/37 j* (f xv^^A^ufw A Cxv^ufw) f w $ 'lpyjdTTl w 
Mzs = /(rJC e A<fv /e + C e $ ufe ) T Y$i>edm e 

M 4 i = Mj T 4 , M 42 = Ml, M 43 = Mj 4 , M 44 = f$Z U ’$™dm w , M 45 = 0 
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M 46 f^uw (r w C w % f u , + C v ,v f V) *f>Ti,f u ,)dni u , , A /47 J^uw^'w ^ ipwdrn w > A /43 0 

A /51 = A/j^ , A /52 = A/^, A /53 = A/^, A /54 = ML M 55 = f $j e $ ue dm e 

M 5& = f$Ie(L C e$t>fe + C e fj e ^, fe )dm e , A/ 57 = 0, A/ 58 = f$l e rj$^ e dm e 
A/61 = A/^, A/62 = A'/Jg, A/63 = A/ : ^ , A/64 = A/^, A/ 6 5 = A/56 

A/66 = ^ipfdnif 4 * f (f u ,Cw^n'fw 4 - Cu' ff hi % fw ) ( f w ^xl'fxv A - C W T j w ^-^f w )d'l7l w 

+ j if e Ce^rpfe 4" ^e^/e^/e) (^e ^e^V/e A" E e f j- e ^^,f e )dTn e 
A/57 = J (r w C w $tf W 4 " Cw'r ^ % v^iirwdTn w 
A/es = j(rjc e % fe + C e fJ e ^^fe) T L $^dm e 

A/71 = AC A/72 = ML A/73 = A/ 3 t 7 , A/74 - Ml, A/75 = A/57, A/ 76 = A/ 6 r 7 
A/77 = //?8 = 0 

A/si = A-Z/g, A/ 8 2 = A/^, A/ 8 3 = A/^ , A/ 84 = A/^, A/ 85 = A/^, A/ 86 = ML A/ 8 7 = MT 8 
Mss = J^l e f e fJ $ n , e dm e 

in which m is the aircraft total mass, 

S — 1 (r j 4" ^ufO.uf'jdTTlj 4” /' [( Ay 1 /' 4~ A u f u:^\uf ) ic 4" C w (f w 4" u , ) 1 lr (i 

+ J[if fe + ^u/eQu/)^ + Cf ( f e 4- <& ue (\ ue )\C e dm e (16) 

is the matrix of first moments of inertia of the deformed aircraft and 

J = jiff A - ^u/Qw/) iff A” *&ufQuf)dmf + f[C w (ff w + ^u/ioQu/) 4* if w $ uu ,cj uu i) Cw] 

X [C w (r f w 4" &ufu’Quf) A - ifw A - ^ UW Quit’ ) Cif^dlTTlyj “I - jiC e (j'^ e “I - *&ufeQuf') 

4- (f e A- $ ue q ue) T C e ] T [C e (ffe 4- $ufe<\uf) T + if e -f $ ue q ue ) T C e }dm e (17) 


is the inertia matrix of the deformed aircraft. 

In a similar fashion, we insert Eqs. (10) into Eqs. (9), integrate over the respective components 
domains and obtain the generalized Rayleigh’s dissipation functions 


„ d 2 $ T 4 2 <b . 

d~ia — Jf) i d~ ut dD, = | f^C-ujEI, q u j ~fff2 d x 2 QuidDi — ^Qm/'tnQui > 

\pi ~ , f[) : f ~ !.’J dD : — - Cyji G ■/, q 2 y ~dx Ol' d Dj = 2 Qi.; • 


i = f,w,e (18) 


where 


C u i j D c U iEf 


d?$l d 2 ^ 


dx 2 

l 


dx 2 


~ dD t , Cjpi J y-j c^jj (? / j 


d&l d%j 

dx, dxi 


dD it i = f,w,e (19) 


are damping matrices. 

Next, we denote the momentum vector for the whole aircraft by p = [p^ / pj f pj f pf w pf e p^ f 

pL PleY = [pf pI • • • pJL S0 that we can write 


p = dT/dV = A/V 


( 20 ) 
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where the individual momenta are given by 


8 

pvf = 9T /av, = p, = Yi m^Vj 

3 = 1 
8 

Pa-/ = dT/du>f = p 2 = J2 M 2 jVj 

3 = 1 
8 

P uf = dT/ds U f = p 3 = ^2 m 3 jV j 

3 = 1 


( 21 ) 


8 

Pv-e = dT/ds^e = p 8 = ^2 M s jVj 

3 = 1 


Finally, adding some obvious kinematical identities to the discretized version of Eqs. (1), the state 
equations can be written in the special form 


R / = CjVf, Of = Ej'utf-, q ui = s ui , q ^ = s^, i = f,w,e 
Pvf = -WfPvf + F, = -VfPvf - W/P w/ + M 


P ui dT/d q ui - 

KuiQui -c ui 

&ui + Q ui 

Pxpi q^'? "f” 


i = f,w,e 


(22) 


where, using Eq. (8) in conjunction with the second form of V/, V^, and V e , Eqs. (12), 

dT dVj dT dVl dT dVT dT 

dq uf dq ufdVf dq uf dV w dq uf dV e 


= f$Z f u;jV f dm f + $Z fw rfCZfV w dm w + fe uJC[/V e dm e 

dT dVj dT dVl dT dVT dT ~ t _ 

« — = + fl m 7 + Q = *„C bW/ jy w dm w 

&Q.UW &Quw oV f UQxlw d\ w dQiiiu 0\ e 


dT 

9q U e 


dVj dT dVl dT dVj dT 
dq ue dV f dq ue dV w dq ue dV e 


= $l e C e Vf T JVedm e 


Moreover, 


E u i — J^uiEui^uidDj, K- l i n — i f, w,e 


(23) 


(24) 


are component stiffness matrices. In practice, K ui and K^i ( i = /, w, e) can be obtained with 
greater ease from the strain energy directly, as shown in the Numerical Example. 

The quantities F, M, Q ui and (i = /, w, e ) appearing in the state equations, Eqs. (22), 
represent generalized forces. They are related to the actual forces, which consist of the distributed 
force fi(rj,<) over component i due to gravity, aerodynamics and controls and the engine thrust 
F E S(r — r E ), where <)(r — r^) is a spatial Dirac delta function, 37 in which r E denotes the location 
of the engines. If some control forces are concentrated, they can also be treated as distributed, as 
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in the case of the engine thrust. The relation between the generalized forces and the actual forces 
can be obtained by means of the virtual work, which can be expressed as 

SW = Y« f Dl \tf + F£S(r - r^dR'dA (25) 

where 5R* is the virtual quasi-displacement vector of a typical point on component i (i = /, w, e). 
The vector &R* is related to the virtual quasi-displacement vectors corresponding to the quasi- 
velocities used to describe the motion of the aircraft components. Indeed, using Eqs. (12), we can 
write 

5R} = 5RJ- + if/ + $> U f(i u f) T 50*j: + ^ u fdci u f + 

<5R*, = C' tl ,5Ry + [C w (ff w + ^ > ufw ( luf) T + (^w + < &uw ( luu’) r Cu<] / 5@f (26) 

+ (r^,C w Ad? u f w + C„,$ u / u ,)(5q u / + $ uw^Quw 

%■ f w ) “b ^ ii’u' ^qi. ■ w 

and we note that 6R* can be obtained from <£R*, by replacing w by e. Inserting Eqs. (26) into Eq. 
(25), and collecting terms, we can write the virtual work in terms of virtual generalized displace- 
ments as follows: 


sw = F t 6K} + M T 5d) + ZMMui + (27) 

from which, assuming that the engines are mounted on the fuselage (Fig. 1), we obtain 

F = L,llj + Fe<5(f - T E ))dD, + CZf D J w dD w + CjS D f,dD e 

M = J D f(rf-{-<t> u fC[ u f){ff +F£;5(r — T£)\dDf + Jp w [(^fw + *&ufwQ.uf )C W 

T" (Tin ~b *&uu’Quw^\^wd'Du! "|- J [(f'/ e “b ^ufeQuf -b C e (t e $ ue q ue )]f e d77 e 

Q u / = J D $l f [if + F £ 5(r - T B ))dDf + J D (f^C u ,A4> u/u , + C w <b ufw ) T f w dD w 

W (28) 
+ f De (r JC e A<f> U f e 4- C e <$> u f e ) T f e dD e , 

Qt/ = ~b F/-;() (r — Y //))(! Df + f £) u {^yCu'^jpf-w "b fj pdD w 

+ f ^(rJCeQ^fe + C e fJ e $Tp fe ) T f e dD e 

Qui = JD®uA dD ii Qipi = Jd^^^Du i = W,e 

To complete the state equations, Eqs. (22), it is necessary to derive the stiffness matrices K ui and 
K^i (i = /, w, e), the aerodynamic forces and the control forces. 
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4. The Aerodynamic and Gravity Forces 


The forces acting on the aircraft can be identified as the aerodynamic, propulsion, control and 
gravity forces. Equations (28) give the generalized forces in terms of actual distributed and con- 
centrated forces, where the first imply prescribing the forces at every point of the aircraft. Of 
the forces acting on an aircraft, the aerodynamic and control forces require special attention. We 
discuss the aerodynamic forces in this section and the control forces in a later section. 

There are a number of aerodynamic theories available, some of them capable of prescribing the 
pressure distribution at every point of the aircraft. However, as pointed out in the Introduction 
and Literature Review, any aerodynamic theory to be used in a dynamic simulation such as that 
described in this paper must satisfy certain requirements. Indeed, one of the requirements is that 
the aerodynamic forces be expressed in a form compatible with the present general dynamic for- 
mulation, which implies that they be in terms of the same variables and be referred to the same 
body axes attached to the undeformed aircraft as here. Another requirement is that the aerody- 
namic forces lend themselves to sufficiently fast computation as to permit time simulation of the 
aircraft behavior. Such computational efficiency does not appear to be within the state of the art. 
Hence, an aerodynamic theory capable of fitting within the framework of a computer simulation 
of the type envisioned here must yet be developed. Such a theory need not be unduly accurate, be- 
cause a feedback control design tends to be sufficiently forgiving to tolerate small deviations from 
the actual aerodynamic forces. Until such a theory becomes reality, it is still possible to demon- 
strate how aerodynamics fits in the integration process by using an existing theory satisfying the 
requirement described above, namely, strip theory, 38 even though the theory may not be entirely 
suitable otherwise. 

To derive the aerodynamic forces included in the distributed forces f (i = f, u>, e) acting on the 
aircraft by means of strip theory, we regard the fuselage, wing and empennage as two-dimensional 
aerodynamic surfaces. The lift force per unit span of fuselage can be written as 38 


if = QfCfCLafOtf = qfCfCLf (29) 

where cj is the chord, C La j the slope of the lift curve, C L f the lift coefficient and 

q f = \p(V% + V f l), a f = tan ~\V f jV fx ) (30) 

are the dynamic pressure and angle of attack, respectively, in which p is the air density and Vj x 
and Vf z are components of the velocity vector V /, Eq. (5). Similarly, the drag force per unit span 
of fuselage is given by 

df = QfCfC D f = qfCf(C Df0 + kfCl f ) = q f c f (C D fo + kfCl af a? f ) (31) 

where Cp/o is the drag coefficient corresponding to aj = 0 and kf is a constant. The fuselage 
has also vertical surfaces subjected to aerodynamic forces. The lateral force per unit span can be 
expressed as 


S/ QsfCsfCspfPf — Qs/CsfCsf (32) 

where c s f is the lateral chord, C s pf the slope of the lateral force curve, C s f the lateral force coeffi- 
cient and 

l.f = + «%), /»/ = tan-'ftyv,*) (33) 
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are the dynamic pressure and the angle of attack of the lateral force, respectively, in which Vj x and 
V fy are components of V /. Hence, the aerodynamic forces on the fuselage can be written in the 
vector form 



i f sin a f — df cos aj 


Sf sin (3f 


0 

? f »/ = 

— Sf cos (3 f 


— t f cos ctf — df sin a f 


0 


(34) 


In a similar fashion, the lift and drag per unit span of wing are given by 


{C LctwQw~^- & LSa^a) i Qw^w Lw) DwO~^~ Law^w ) (35) 


where c w is the chord, S a an aileron rotation, C L s a a control effectiveness coefficient 1 ' and 


Qw = \p(Vwy + Vwz)i a w = tan l (V wz /V wy ) 4- Vw (36) 

in which the velocity components V wy and V wz of V w can be obtained from Eq. (6). Moreover, 
ip wx is the angular displacement of the wing about axis x w due to torsion. There is no meaningful 
lateral force on the wing, so that 


0 

f aw | sin OCyj dw cos o „ 

—t w cos a w — d w sin a v 


(37) 


The empennage has both lift and lateral surfaces. The lift, drag and lateral forces per unit span of 
empennage are 

-{■" C LSe^e) ^ dg Qe^ei,^ DeO ""f ^ e L&c ) 1 QseCsei^* sfiefte "i - ^ sSr^r^) (38) 


where c e and c se are the chords, 5 e and S r are rotations of the elevator and ruder. Cue and C st s r are 
respective control effectiveness coefficients and 

Qe — \p(Yey V e z)i & e = t an {Vez/Vey) + V* ex ( 39 ) 

Qse = W? y + V e l), P e = tB.Xl-\V ez /V ey ) + </>e* 

Hence, the aerodynamic force vectors per unit span of empennage can be written as 



0 


0 

f = 

L ae 

sin (fg cos ftg 

f = 

) L se 

5 e sin fie 


^g cos Oig c(g sin dg 


_ S e COS fie _ 


(40) 


For a typical component, the lift, lateral force and drag per unit span are applied at the line 
of aerodynamic centers. Hence, in Eqs. (28), the domain of integration for the terms involving 
the aerodynamic forces is the line of aerodynamic centers. The gravity forces per unit volume of 
components are simply 


0 


0 


0 

0 

? ^gw — CwCf 

0 

, fge = CeCf 

0 

. Pf9 . 


PwQ 


. Pe9 . 


( 41 ) 


Note that the aerodynamic and gravity forces are in terms of respective component body axes. 
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5. A Unified Theory for Flight Dynamics and Aeroelasticity 


To complete the discussion of the forces acting on an aircraft, we turn our attention to the control 
forces. Aircraft control is carried out by means of control surfaces, as well as by the engine thrust. 
Before we consider the problem of control design, it will prove beneficial to examine the nature 
of the controls. Controls are of two general types, one designed to steer the aircraft as a whole 
and the other to suppress the effects of any undesirable disturbances. The first type involves rigid 
body motions of the aircraft, which are in general large, and traditionally lies in the domain of 
flight dynamics. On the other hand, the second type involves elastic deformations of the aircraft, 
which tend to be small compared to the rigid body motions, and traditionally lies in the domain of 
aeroelasticity. Hence, the formulation given by Eqs. (22) can be regarded as spanning the fields of 
flight dynamics and aeroelasticity. 

From the above discussion, it appears that control of the aircraft as a whole is likely to be 
different in nature from control of disturbances. In this regard, we observe that the state equations, 
Eqs. (22), are in general nonlinear and of high dimension, where the nonlinearity can be traced to 
the large rigid body variables. On the other hand, the high dimensionality can be traced to the small 
elastic variables. In view of this, a solution by a perturbation approach seems indicated, which 
amounts to a separation of the problem into a zero-order problem for the large variables and a first- 
order problem for the small variables, where the difference between the large and small variables 
is one order of magnitude, or more. Physically, in the zero-order problem the aircraft executes 
a given maneuver as if it were rigid, in which case the mathematical formulation consists of a 
maximum of six coupled, generally nonlinear ordinary differential equations, three for rigid body 
translations and three for rigid body rotations. They correspond to the equations commonly used 
in flight dynamics. On the other hand, the first-order problem involves the elastic deformations, 
as well as small perturbations in the rigid body variables. In view of the inclusion of the latter, 
the first-order problem defined here represents an extended aeroelasticity theory, in which the 
rigid body degrees of freedom are included in a natural manner. This is in contrast with some 
occasional practice, 11,39 in which “body freedoms” are included in an ad hoc manner. We note that 
the solution of the zero-order problem enters into the first-order problem, so that this new theory 
provides one set of extended aeroelasticity equations for every conceivable rigid body maneuver 
of the aircraft, rather than the single set of equations commonly associated with steady cruise. We 
express the perturbation solution as follows: 


R f = R<°) + Rf, e, = ef + ef. v, = vf + vf, u , = <*f + Jf 

Pvt = Pvt + Pv/. Pvt = pi°) + pi/, F = F<°> + F<‘>, M = M<°> + M'» 

where the superscripts (0) and (1) denote orders of magnitude. All the quantities related to the 
elastic deformations are regarded as being of first order. Then, inserting Eqs. (42) into the state 
equations, Eqs. (22), and separating different orders of magnitude, we obtain the zero-order prob- 
lem, or the flight dynamics problem 

Rf = cf T vf , ef = 

pf, = - Zf pj?> + F<°>, pf = -Vf >pf - of pf + M<°> 

in which Cj 0) and E { f 0) can be obtained from Cf and Ef, Eqs. (2), by replacing ip, 9 and 4> by 

°), 9 (0 ^ and f (0) , respectively. Moreover, from Eqs. (28), the zero-order generalized force and 
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moment are given by 

F (0) = Lj|f (0) + F m J(r _ TE)]dDj + clS D f»dD u , + CJS D t^dD. 

M<°> = f D /,{tP + F fS(r - r E )] dD, + f D Jf + Clr^dD, (44) 

+ f D S?H c T + cT?')tl 0)dD e 

where the zero-order parts of the aerodynamic force densities contributing to f/' 1 , fu 0i and fi 0) are 


so 

II 

iP sin — dp cos q| 0) 

0 

—ip cos ap } — dP sin ap 

II 

o' ^ 

'"t" 

tM 

sP sin pp 
—sP cos pp 

0 


f(0) _ 

l ai 

0 

ip sin — dP cos a) 0 ^ 

—iP cos a| 0) — dP sin a| 0) 

, i = w,e; fP = 

si 01 

L -4 o) 

0 

1 sin dP 
cos dP 


in which the zero-order parts of the lift, drag and lateral forces per unit area are 

tf=qfc f C Laf a { p, C=q { ^c w {C Law a^+C L6a 5^), 4°> = qPc e (C Lae aP + C LSe SP) 
sf> =q < ffc sf C s0f (3f\ 4°) = qi?c se (C s0e pW + C s5r 6i 0) ) (46) 

df ] =qf ) c i [C Di o + kiCl ai {a \ 0) ) 2 ], i = f,w,e 


where 

q?> = 'iP{(V}y + (Vj!”') 2 ], a<°> = tan-'(V';“ , /F / ( ">), /3<°» = tan -'(Vff/V™) 

' =ip[(V 4 J, 0) ) 2 +(V;t 0 >) 2 ], Qf^tan-'IV'f/V’®). /f ’“tan- 1 ^ 01 /^’). i = W?) 

$ =§p[(f / , ° > ) 2 + ( v^n = jphk ', 01 ) 2 + (h™ 1 ) 2 ] 

Moreover, the zero-order parts of the gravity force densities are 


C = cf [o 0 P/ p] r ; = acf [0 0 Pi <?] r , i = w, e (48) 

The zero-order state is defined as x (0) = [R^ 0)r 0 ( ® )T pPp ppp) T , and Eqs. (43) contain in 

addition Vp and Pp. However, V^ 0) and u>f ] can be expresed in terms of p[?j and pP f by using 
Eqs. (21) and writing 

pW = m vf + S (0)T Jj ?\ p $ = S {0) Vf ] + J<°) uf (49) 

The solution of Eqs. (43) in conjunction with Eqs. (44)-(49) consists of the state x (0) and it repre- 
sents a given maneuver of the aircraft. 

Inserting Eqs. (42) into Eqs. (22) and retaining the first-order terms, we obtain the extended 
aeroelasticity problem defined by 


it' 1 ’ = cfvf + c} ,,T vf , ef = ( Ef y'jp - 

() u x — &ui ■ Q ipi ~ S VJ > i f ,W, € 
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*(1) _ _ ,~,(1)„(0) _ .flJ 1 ) , p.(l) 

P V f — ^ f P Vf ^ f P Vf ' 


p s - - vj'V* - if pg> - - fi>J»pS» + M»> 

P»/ = /^/(^V/ 11 +«} 1)r vf Mm/ + K,J“f >T CZVl ;> 

+ f*l„(>f >T CW + iif T CjVf')dm e - K„ q„/ - C„,s u/ + Q u , 

+ / < I > I/i/ 0,T V™Mm / + f$l f ,^ , ° lT ClV'"'dm u , + f$l le uf ,J Cj V“”dm, 
P t/>/ — “ ~ ^W S ^/ + Q^/ 


P«f = /®»i(Ci"/ 0> VjP+Ciw ” 1 Vf ) )rfm i -A'„ i q„ i -C„ i s„ i +Q„ j 
+/^C ( u.f T V< 0 , <Jm j , »=«o,e 

PV'i ^ ipi (-'ijji&ipi Qxpi i 2 — Ti) ? C 

where 

c} 1 * - c}°V (1) + cj°> (1) + c^ {1) , £} 1} = Efje {1) + £<°V (1) 

in which 


(o) _ dC f 


dip 


^ ( 0 ) ,^ 0 ) ,<a (0) 

— sip^cO^ cip^cO^ 0 

— s^ 0 )s 6 ^ 0 )s<^ 0 ) — cip^ c<j>^ cip^sQ^scp^ — 0 

—sip^sd^c<f>^ + cip^scp^ -f s'tp^scp^ 0 


r (0) _ 9C, 

^ fa — 


'19 


MO) _ 

^ “ d<t> 


( 50 ) 


(51) 


d® lv>(°),e(°), 0 <o) 

—cip^s 0 ^ —sip^sd^ — c 0 ^ 
cip^ c 6 ^ s'tp^ cd^ s<p^ —s 0 ^s <fi^ 

CW i °' > c9 { - O ' ) C0 (O) Slp^o 9^C(f)^ —S0^C :0(°1 

8C/I 

li/>< o >, 0 < o >, 0 < o > 

0 0 0 

S'(jj^s9^C <f>^ — s4>(°^ C0^C </>^ 

—Clp^sO^sfi^+SIp^C^ 0 } —S1p^S0^S(f)^ — CXp(°} C<j>(°) —cd^Stff^ 0 ) 


(52) 


r(°) - 
J f 0 - 


dE f 


39 


e<°) , 0 ( 0 ) 


(0) _ dEf 
^ - d<j> 


0(0) ,0(0) 


0 0 -c 0 ( o) 

0 0 — 

0 0 -s0 (o) c0 (o) 

0 0 0 

0 —s ^ 0 ' 1 cd^c<p^ 
0 —c(fi^ —c 9 ^s <j>^ 
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and 


V / 0) = V / 0> + V/* = Vj: 1 * + rjo;^ + $ u/ q u/ r u>^ 0) + $ u/ s u/ + rj^fS^f 

V< 0) = C u ,Vf + (0*7% + 

v£> = c w vp + (c w fj w + fZc w )uf + (c u ,$ t( ^q u / + ( 

H - )^w/ “t” 

“I - ( C'W f w ^Xpfw )®V ; / "f" ^*U» ^ *l' w 

Then, from Eqs. (28), the first-order generalized forces are 

F (,) = Id, l ff +F^(r - r e ))dD, + CZJ D /J>dD„ + Cj f^dD, 

M (1 > = + F™<5(r - r E )] + + Fg’f(r - r E )]}<i£>, + f D Jf /w Cl 

+ Cjr.)4‘» + + C'j4.^,„„)f™»]<iA,. + i„,[(r/ t C-J + Cff e )f<» 

+ (<^q k/Cf + 


Q«/ = / D ,C/[f/ 0> + f/ 11 + Fi?<5(r - te) + Ff6(r - T E )]dD, + J D JflC w At„ w (54) 

+ c w $ ufw ) T (C> + 1 l»)dD w + f D jf!c,A* u „ + c e $„„) T ( f<°> + r:">dD, 

Q*1 = f D *1, fll^ + + Fg>f(r - r E ) + F»><5(r - r E )]dD, + / D .(fJC»4. w „ 

+ C w f] w <S> t , w ) T ( tS>+tS')dD. + f D '(flC'i He +C e fT e % fe ) T (t<»+y>)dD, 

Q«, = f D| t>Zi(f! 0, + <fV*. Q.« = + fl'IdD,, i = w,e 


where the first-order parts of the aerodynamic force densities contributing to f^ 11 , fi 1 * and f,: are 


r(l) 






Pp sin Qj 0) — d!p cos ap + pp cos q| 0) 

0 

— ^ 1) cosa| 0) — c?^ sinay 0) + pp sina^ 

s { p sin f3 ( p + Sf ] f3 { p cos (3p 

— cos f3p' 1 + flP sin / 3 P 

0 


+ d,p sin *P)af 
— dP cos ap)aP 


0 

#P sin ap — d ( p cos af^ + (/p cos ap 
—d-P cosa- 0) — dp sina) 0) + {iP sin a- 0) 


+ dP sinaf^aj^ 
— dP cosa| 0) )a)^ 


0 


sP sin ftP + sPdp cos dP > 
sP cos ( 3 p + sPfdP sin i 3 i 0> 


(55) 


i = w,e 
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in which the first-order parts of the lift, drag and lateral forces per unit area are 

cj — CfL>i a f(qj Ckjr + Qj Otf ) 

4 X) = q^UCLa^ + CuaSP) + i° ) C w (C Iaw ag ) + CuJP) 
4 X) = q^c e {C Lae a^ + C LSe Si 0) ) + qi 0) c e (C Lae a^ + C LSe S^) 
s f ] = c sfC S 0f{q^jpf ) + qffpf) 

4 } = q^CseiC^ + C sSr S< 0) ) + q ( s °Jc se (C s0e p^ + C sSr 6W) 


( 56 ) 


rf! 1} = ql 1} Ci[C Di o + k r Cl ai (af ) ) 2 } + 2qf ) c i k i Cl <xi af ) af\ i = f,w,e 


where 


,<*) 


qy = p(v}°>vff + Vj"’Vjl>), af = tan" 1 

= Mvffvg’ + V^iy), 0? - 


vjyff-vgvn' 


(v f fr + (v f '"y 

/a: v fy 


K0)> 


tan 


-l 


yWym _ fWy<oy 


(V<y + (v;;<) 

y-Wy-U) _ y(l)y(0) ' 

21/ 22 r Zj/ *22 

(f^) 2 + «i 0 ’) 2 


(°)\2 


q? ) =p{vS ) v" + V l 


(0) V^ 1) ), aj x) = tan" 1 


(57) 


«2> = rtKSW + = tan- 1 


y(°)y (1) _ y 
^et/ *e2 re 


(l)f>(0) 
re.z 


ej/ 


+ ipix, i = w,e 


+ tpe 


(Ve?) 2 + (Ve^) 2 

Moreover, the first-order parts of the gravity force densities are 

$ = Cf [0 0 p f g ] T ; #> = CiCf [0 0 Pi g) T ,i = w,e 

Equations (50)-(58) must be solved for the state x (1) = [K^ )T 0^ T q^ f q^ u , . . . q£ e 
P 1/ PL • • • pJeF in conjunction with 


(58) 

(1)T (1 )T 

P Vf P If 


pi 1 } = mVf + S (1)r u>} 0) + 5 (0)T ^} 1} + M^suf + M'S s uw + M$ s ue + Mjfs*/ 

+ M[jSip w + gSfe 

p(l) = §(l)yW + 5(0) V (1) + J(1) W (0) + J(0) W (D + M (0) Su/ + M (0) Suw + M (0) 

+ MS Sip f + MS Syw + M^Sxpe 
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P., = Mg’v<” + Mg’V™ + Mg’u,}” + Mg’w™ + MS’s., + MS’s... + Mg 

“I" ^3 6 s ipf + 

P™ = + M^vf + M^u,^ + s ue 

■I" + M^Sipe 

p„e = Mg’V* 1 ’ + M<;>V™ + Mg’wj,” + Mg’wf + Mg’s., + Mg’s.. + Mg’s« 
+ Mg’s#, + Mg’s*. + Mg’s*. 


r(i)v(°) 


r(0), .(1) 


(1), ,(0) 


r(«». 


r<°). 


r(°). 


(59) 
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P */ = MgW 1 " + MgWf 1 + M^uf + Mg' uf + Mg' s„ + Mg's m + Mg', 

, 71 ,r(0) , n /rfO) , n ^(0) 


*ue 


+ + Mffi) Sip w + Mgg Sip e 


P*w = + M 7 ( ; } Vf + + Mg ] ujf + M^s u} + M^s^, + Afffs, 

'(®)« I 71/f( 0 ) 


A/ig } s^ + Mf^s 


_!_ TLfW 
'tyw I -‘■VJ-'fg 


p^ e = M&tyV + M£>Vf ] + M<° } uf + M$? 

_l 


+ Mor s^f + Mo-) i 


iS°) I I Jlf^a i 

U7 f + M 83 S U J- + iW 84 s uu . -+- m 85 s, 


where Af^ 0) and are the zero-order part and first-order part of M tJ , Eqs. (15). Equations (59) 
are to be solved for VjV, u> { j\ s uf , s uw , ... , s ve in terms of p- 1 }, p ( J}, p uf , p uw , ... , p te , Vj 0) 
and and the result inserted in Eqs. (50). 

The zero-order problem, or flight dynamics problem, represents an inverse problem, which 
amounts to determining the controls permitting realization of a given rigid body maneuver. The 
first-order equations representing the extended aeroelasticity problem are linear and tend to be of 
high order. Moreover, they contain the zero-order variables and u^ 0 *, representing a given 
maneuver, as coefficients and as an input. If V^ 0) and are constant, then the system is time- 

invariant, and if V)' 1 and lo)' depend on time, then the system is time-varying. In either case, 
controls can be designed by various methods. In the time-invariant case, a stability analysis for the 
closed-loop system can be carried out by solving an eigenvalue problem. Such a stability analysis 
is precluded in the time-varying case. Simulation of the response of the closed-loop system to 
external excitations, such as gusts, can be obtained in both the time-invariant case and time- varying 
case. 
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6. Control Design 


Flying aircraft are subjected to various disturbances tending to drive them from the intended ma- 
neuver and to cause vibration. If the system is controllable, these effects can be suppressed through 
controls, which are carried out by means of actuators; in the case of aircraft they consist of the en- 
gine thrust and the control surfaces. A system is said to be controllable if there exists a piecewise 
continuous input that will drive the initial state to any final state within a finite time interval. Sys- 
tem controllability can be determined by checking the rank of a so-called controllability matrix, 40 
which may not be feasible if the system order is large. In practice, controllability can be ascer- 
tained on physical grounds by making sure that the input forces, namely, the forces due to the 
engine thrust and control surfaces, affect all the state variables. 

As indicated in earlier sections, there are two types of controls, one type designed to permit the 
aircraft to execute a desired maneuver as if it were rigid and the other type to reduce any deviations 
from the rigid body maneuver to zero, which amounts to suppressing vibration and perturbations 
in the rigid body motions of the aircraft. The first is associated with the flight dynamics problem 
and the second with the extended aeroelasticity problem. In general, the engine thrust and control 
surfaces are designed so as to ensure that the aircraft is able to carry out the required maneuvers, 
as well as to suppress any undesirable disturbances, thus addressing the needs of both the flight 
dynamics problem and extended aeroelasticity problem. 

Using Eqs. (43), we write the flight dynamics problem, or the zero-order problem, in the com- 
pact state form 

x ,0 »(() = f[x |0 >((), VS>(<)] + B (0, [V™(t)] U <0 >(f) (60) 

which, from Eqs. (49), must be considered in conjunction with 

PrtW = Af*v!?(t) <6I) 


where the zero-order quantities are identified as follows: x (0) = [R^ 0)T ppj pi°/ T ] r is 

the state vector, f is a nonlinear function of the state vector and the zero-order rigid body velocity 

vector = [V^p u>P ] T , B ^ is a coefficient matrix, u ( °) = [F^ 0) SP <5e 0) $ 0) ] T is the 
control vector, in which Fp, SP, dP and SP are the engine thrust and control surfaces angles, 
= [pPj pp f ] r is the rigid body momentum vector and 


MP = 


ml SM T ‘ 
5(0 j{ o) 


(62) 


is the rigid-body mass matrix. In the context of the present integrated approach, Eqs. (60) and (61 ) 
represent an inverse problem, in the sense that a state vector x f0) describing a desired maneuver is 
postulated and a force vector u (0) permitting a realization of the given maneuver is determined. 

Next, we assume that x (0) (f) and V^(t) are known and use Eqs. (50) and (59) to express the 
extended aeroelasticity problem, or first-order problem, in the form 

x (1) (f) = A(f)x (1) (f) + B(t)u w (t) + [0 7] r F ext (f) (63) 


and 


p (1) (f) = M (0) V (1) (f) + M (1) (f)V (0) (i) (64) 
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respectively, in which x^ = [R (1)r 


aW T T n T T „(! )T _.(! )T T T „T 1 T 

/ u f Huf VI uw • ' ' t*Vf ”u>f P uf Puw ’ * * P^eJ 

is the first-order state vector, A(t ) = y4[x^°^(t), Vj.^(£)] and £?(£) = B [x^(i), V^(t)] are 

coefficient matrices, u^(t) = [F ( ^ } 6a } 6 ( f l) <5r^] r is a first-order control vector, F ext is an 

external disturbing force vector, such as due to gusts, p (1) = [p^j. p ( JJ p^f p^ w ■ ■ • pjj r is the 
first-order momentum vector and 
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are zero-order and first-order extended mass matrices. Note that F ext is regarded as a disturbing 
force of a transient nature whose effects will be eventually suppressed by the controls. 

Equation (63) represents a set of linear equations, and the objective is to find a control vector 
u (1) (f) that drives the state vector x (1) to zero. To this end, we consider a linear regulator whereby 
the control vector is a linear function of the state vector. In particular, we consider a linear quadratic 
regulator (LQR) in which the objective is to determine an optimal control vector minimizing the 
quadratic performance measure 40 


J = ix (1)r (f/)ifx (1) (t/) -I- \ [x (1)r (f)Q(f)x (1) (£) -I- u (1)r (t)R(t)u (1> (f)]dt (66) 

where H and Q are real symmetric positive semidefinite matrices R is a real symmetric positive 
definite matrix, t 0 is the initial time (commonly assumed to be zero) and tf is the final time. It is 
shown in Ref. 40 that the optimal feedback control vector is given by 


u w (t) = -R~ 1 (t)B T (t)K(t)x {1) (t) = -G(t)x (1) (t) (67) 


where 


G{t ) = R~\t)B T (t)K{t) (68) 

is a control gain matrix, in which K(t ) is a real symmetric matrix satisfying the transient matrix 
Riccati equation 

K = -Q - A t K - KA + KBR~ 1 B t K, K(t f ) = H(t f ) = H (69) 

a nonlinear equation which must be integrated backward in time from tf to t 0 . Rather than inte- 
grating a nonlinear matrix equation, it is advisable to transform the problem into a linear one. To 
this end, we consider the transformation 

K(t) = E{t)F~\t) (70) 

where E(t) and F(t) can be obtained by solving the linear equation 
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-A T (t) 
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E(t,) } 


! H ' 

F(t) _ 
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which again must be integrated backward in time. Inserting Eq. (67) into Eq. (63), we obtain the 
closed-loop equation 

*“> (<) = [>!(<)- B(<)G(<)]* (1) (<) + [0 /] r F„(i) (72) 
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which can be integrated to simulate the system response. 

The problem is considerably simpler when the zero-order solution is constant, such as in steady 
level flight, as in this case the coefficient matrix A is constant. Then, if the system is controllable, 
H = 0 and Q and R are constant, the Riccati matrix K approaches a constant value as t; increases 
without bounds. In this case, Eq. (69) reduces to the steady-state matrix Riccati equation 

-Q - A t K -KA + KBR~ 1 B t K = 0 (73) 

a nonlinear algebraic matrix equation, which can be solved by means of Potter's algorithm, 40 and 
the gain matrix G becomes constant. The closed-loop equation reduces to one with constant coef- 
ficients, or 


x (1) (f) = (A - BG)x w (t) + [0 I] T F ext (t) (74) 

which can be used for response simulation. For a stability analysis, we solve the associated eigen- 
value problem 


(. A-BG - A I)x = 0 (75) 

The closed-loop system is stable if all the eigenvalues are pure imaginary and/or complex with 
negative real part. 

The control vector u (1) is optimal in the sense that it minimizes the performance index J, but the 
physical merit of this optimality is debatable. In fact, it is often necessary to adjust the otherwise 
arbitrary weighting matrices Q and R to achieve a desirable system performance. The real value 
of the LQR algorithm is that it guarantees a stable closed-loop system. 
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7. Optimal State Observer 


In Sec. 6, we carried out the control design in two stages. In the first stage, we postulated a 
desired aircraft maneuver x (0) and used an inverse approach to determine the control vector u (0) 
permitting realization of the maneuver. In the second stage, we designed a feedback control vector 
u (1 ) ensuring stability of the maneuver, which amounts to driving the perturbation vector x (1 ) to 
zero. The control vector is given by Eq. (67), in which G is the control gain matrix. 

Implementation of the control law, Eq. (67), requires knowledge of the state vector x^\ which 
can be obtained through measurement. This creates somewhat of a problem, as measurements 
represent real quantities and our state vector consists of abstract generalized coordinates rather 
than real coordinates. Moreover, we feed back only perturbations from the maneuver variables. In 
the absence of external forces, the state equations describing the extended aeroelasticity problem 
have the vector form 


x (1) (f) = Ax (1) (i) + £?u (1) (t) (76) 

where the coefficient matrices A and B were defined in Sec. 6. Moreover, denoting the measure- 
ment vector by y (t), we can write 

y(t) = y {0} (t) + y^(t) (77) 

where y (0) (f) is the contribution from the maneuver variables and y (1) (f) is the contribution from 
the perturbations in the maneuver variables. We express the latter in the form 

y (1) (t) = Cx (1) (t) (78) 

and refer to y^(t) as the output vector. The assumption is made here that the system is observ- 
able, 40 which implies that the initial state x (1 ^(0) can be deduced from the outputs within a finite 
time period. Observability can be established by checking the rank of a so-called observability 
matrix, 40 a matrix involving the matrices A and C. For large-order systems, working with the ob- 
servability matrix may not be feasible, and in practice the choice of sensors ensuring observability 
must be made on physical grounds. In particular, the choice must be such that the sensors signals 
permit reconstruction of the state at all times. The task of determining the matrix C relating the 
output vector to the state vector is discussed later in this section. 

In reality, the state vector cannot be determined exactly from the output vector and must be 
estimated. A device permitting an estimate of the state vector is known as an observei 40 and can 
be expressed in the form 

x (1) (f) = Ax (1) (f) + I?u (1) (f) + K 0 [y {1) (t) — Cx (1) (t)] (79) 

where K 0 is an observer gain matrix. Subtracting Eq. (79) from Eq. (76), we obtain 

e(t) - [A - K 0 C]e(t) (80) 

in which 

e(t) — x (1) (f) - x (1) (f) (81) 

represents the observer error vector. The objective is to find a matrix K 0 such that the vector e(t) 
approaches zero as t increases. In this case, x(t) x(t) with time. For a time-invariant system. 
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this amounts to determining K a so that all the eigenvalues of the matrix A — K a C , known as the 
observer poles, lie in the left half of the complex plane. In implementing feedback controls, we 
must use the estimated state x (1) (f), because the actual state x (1 ^(f) is not available. Hence, Eq. 
(67) must be replaced by 

u m (t) = -Gx (1) (<) (82) 


One question of interest is how the choice of observer poles affect the choice of controller poles. 
To answer this question, we use Eqs. (81) and (82) and rewrite Eq. (76) in the form 


x (1) (f) = [A - BG}x {1 \t) + BGe(t) 


Equations (80) and (83) can be combined into 



‘ xW(f) ■ 


‘ A- BG BG 

' x(0(f) ' 

1 

e(i) 


0 A- K 0 C 

. e(f) 


(83) 


(84) 


Because the coefficient matrix in Eq. (84) is block-triangular, the poles of the combined system 
consist of the sum of the poles of A — BG and the poles of A — K 0 C, so that the observer poles 
can be chosen independently of the controller poles. 

An optimal observer gain matrix can be obtained by adopting a stochastic approach, leading to 
the so-called Kalman-Bucy filter, in contrast to a deterministic observer known as a Luenberger 
observer. To this end, we rewrite Eqs. (76) and (78) in the form 


x (1) (t) = Ax.™ (t) + £u {1) (<) + v(f) (85) 

and 


y( 1 \t) = Cx w (t) + w(t) (86) 

where v(f) is known as the state excitation noise and w (t) as the observation noise, or sensor noise. 
It is customary to assume that v(t) and w (t) are white noise processes, with intensities V(t) and 
W(t), respectively, so that the correlation matrices have the form 

£ , {v(fi)v T (f 2 )} = V(ti)S(t 2 - <i), £’{w(f 1 )w T (t 2 )} = - U) (87) 

and that they are uncorrelated, so that 

£{v(ti)w T (f 2 )} = E{w(ti)v T (t 2 )} = 0 (88) 


in which E{-} denotes the expected value. 

The stochastic observer has the same form as that given by Eq. (79) in which the optimal gain 
matrix is determined by minimizing the performance measure 

J 0 = E{e T (t)U(t)e(t)} (89) 

where U ( t ) is a symmetric positive definite weighting matrix. From Ref. 40, the optimal observer 
gain matrix is given by 

K* 0 {t) = Q{t)C T W-\t) (90) 

where Q(t) is the variance matrix of e(t) satisfying the transient matrix Riccati equation 

Q(t) = AQ(t ) + Q(t)A T + V(t) - Q{t)C T W-\t)CQ(t), Q( 0) = Q 0 (91) 
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In the time-invariant case, Eq. (91) reduces to an algebraic matrix Riccati equation yielding a 
steady-state optimal observer gain matrix. The main problem in implementing a Kalman-Bucy 
filter lies in the selection of the noise intensities V(t ) and W(t). 

At this point, we turn our attention to the determination of the matrix C defined by Eq. (78). 
To this end, we must first specify the measurement vector y^(f), which in turn depends on the 
sensors used. In inertial navigation 41 , which represents the process of determining the position 
and attitude of a moving vehicle from self-contained inertial measurements made on board of the 
vehicle, the system consists of a platform containing accelerometers sensing translational motions, 
gyroscopes sensing angular motions and a computer capable of integrating the sensors signals 
to generate the state. Inertial navigation is widely used for aircraft, in which case there are two 
accelerometers aligned with the North and East directions and three rate gyroscopes with the spin 
axes aligned with the North, East and zenith directions. To ensure that gravity does not contaminate 
the accelerometers signals, the platform is made to rotate continuously so as to remain normal to 
the local vertical. The vertical position of the aircraft is measured by means of an altimeter. In 
view of this, we can assume that the system measures the vectors Ry and Of giving the position 
and attitude relative to axes XY Z. In this regard, it should be mentioned that we referred to axes 
XY Z as inertial, but in reality they represent earth-fixed axes. If the current formulation is used to 
describe relatively long flights, then proper allowance must be made for the rotation of the earth. 
Because the same process can be used to determine the zero-order position vectors R (0 ) and 0^°\ 
we can assume that the measurement system is capable of yielding R^ and 0 { p . 

The above process can be used to measure the perturbations in the rigid-body motions of the 
aircraft and the question remains as to how to measure the balance of the variables, namely, the 
elastic variables. To this end, we assume that there are N t (i = /, w, e) sensors measuring veloci- 
ties at the points P k (k = 1,2, , Ni ) of the aircraft components and express the output vector in 
the form 


yl»(t) = [R< 1)T ef T 0 T 0 3 '...0T+[y« 1 >’(«) yg^t) yg>W 


where 


vr(p 2 ,i) ... vr(p„,i)F, 


r(l )T, 


r(l)T 


(1 W( 


(92) 


(93) 


in which V^(Pfc, t) are vectors of velocity measurements. From Eqs. (12) and (53), we can write 

V< 1} (P fc , t) = Vf\t) + f T f {P k )J}\t) + Cb ( p{t)$ uf {P k )ci U f{t) 

+ $u/(Pfc)s u /(f) + fj(P k )% f (P k )s,pf(t) 
d (1) (f) 


= block-diag Cj k 


V^{t) 


, * = 1,2,... ,N f 


(94) 


where 


block-diag C/ k = block-diag[0 0 u )j°^4> u /(Pfc) 0 0 0 0 0 

I fj(P k ) $ uf (P k ) 0 0 f T f{P k )%f{P k ) 0 0] (95) 

and d (1 ) = [R { f 1)T 0 { f 1)T <£ f qL q£ qj/ v<1) = [Vy )T uf T s l f s r uw s[ e sj, 
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s i> w s te\ T - Also fr° m Eqs. (12) and (53), we can write 


v£\Pk, t) - C w \f + (( C w f T {w + fl(P k )C w )Jp + c w u>f ] $ ufw q u/ 

+ CwOJ f ^uw{Pk)f\uw ” 1 " {j'uXP^)^w^^ufw C w ^ufw)^uf 

H - ^uw(Pk)Suw “I" (j* w (Pk)C w & i>fw + C w r f w *&ii>fw) s i>f 

+ fZ(Pk)$ipw{Pk) S ipw 

= block-diag C wk y(i/^ > k = 1,2,... ,N W (96) 

where 

block-diag C wk = block-diag[0 0 C w z/-p$ u f w 0 0 0 0 

^*nt> (-P fe) 0 ( ^ ( Pk ) *&ipfw + C w 

rl(Pk)*+.(Pk) 0] (97) 

In a similar fashion, we can write 

VM(P k , t ) = block-diag C ek , k = 1,2, . . . ,N e (98) 

where 

block-diag C efc = [0 0 C e uj ( p<& u f e 0 C e ujf } $ ue (P k ) 0 0 0 

C, (C e r% + fl(P k )C e ) (fl(P t )C.A<b„ te + C e $ U f.) 0 

(rf (P„)C.% U + 0 ff (P k )^(P t )] (99) 

Equations (94), (96) and (98) are in terms of d ( 0 and V ( 0. They can be transformed into 
expressions in terms of the state = [d^ T p^ T ] r by considering Eq. (64). To this end, we 
observe that M {1) = M (1) ( q u/ , q uw , q ue ), so that we can write 

M(i) V ( o) = My d ( i) (100) 


where Mv is a matrix depending on V <0) . Hence, using Eq. (64), we obtain 

V (1) = (M { 0 ) ) - 1 (p (1) - M v d (1) ) 

so that 


' d (1 ) 1 

I 

0 

y(0 “ 


(M (0) ) -1 _ 


(101) 

(102) 
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Finally, using Eqs. (92), (93), (94), (96), (98) and (102), we conclude that 


'I 0 0... 0 

0 10 ... 0 

block-diag C /i 
block-diag C /2 


block-diag C jn s 
block-diag C w \ 
block-diag C w 2 


I 0 

_(jVf(°))- 1 M v (M (0) ) _1 


block-diag C w n w 
block-diag C e i 
block-diag C e 2 


block-diag C e N e 


(103) 
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8. Numerical Example 


The flight of a flexible aircraft is fully described by Eqs. (43), (49), (50) and (59). A solution of 
these equations requires the aircraft geometry, the mass and stiffness distributions and the aerody- 
namic coefficients. Information pertaining to an actual aircraft was made available by an aircraft 
manufacturer and is listed in the Appendix. 

The solution of the flight dynamics equations, Eqs. (43) and (49), requires the matrices C|° , 
E { p , the total aircraft mass to, the matrix S (0 ^ of the first moments of inertia of the undeformed 
aircraft and the inertia matrix J (0) of the undeformed aircraft. The matrices C ( p and E { p can 
be obtained from Cj and Ef, Eqs. (2), by simply replacing rj>, 6 and 4> by xp(°\ and <^ 0 *, 
respectively. The aerodynamic forces are given by Eqs. (45) and the required coefficients are given 
in the Appendix. Other data required is as follows: engines locations, r^x = [—108.62 37.0 
— 13.96] r in, r # 2 = [—108.62 — 37.0 — 13.96] r in; total aircraft mass, to = 33.5896 lb • s 2 /in; 
matrix of first moments of inertia and inertia matrix 

0 -134.6827 0 ' 

134.6827 0 0 lb • s 2 

0 0 0 

J (104) 

183183.4257 4.7745 -37624.9453 

566328.8970 81.1583 lb • in • s 2 

symm 704218.5794 

Some of the above quantities involve the matrices of direction cosines between the various com- 
ponents body axes and the fuselage body axes. These matrices are listed in the Appendix. The 
flight dynamics problem essentially consists of setting the control surfaces and the engine thrust 
for a given aircraft maneuver, of which the steady cruise is a special case. It essentially amounts 
to the problem of “trimming” the aircraft. As far as this paper is concerned, it provides the input 
V^ 0> , u p ( yj and p^J to the extended aeroelasticity problem. 

The solution of the extended aeroelasticity 
equations, Eqs. (50) and (59), requires an explicit 
choice of the structural model for the aircraft of 
Fig. 1. As a first approximation, the fuselage, 
wing, and both the horizontal and vertical sta- 
bilizers in the empennage are modeled as beams 
clamped at the origin of the respective body axes 
and undergoing bending and torsion. The fuse- 
lage undergoes the bending displacements Uf y 
and uj z and the torsional displacement tpf x , as 
shown in Fig. 2, so that u/ = [0 uj y Uf z ] T and 
\j>f = [ip/x 0 OF. On the other hand, the wing 
and the stabilizers undergo only one bending and 
Figure 2. Aircraft Components Undergoing one torsional displacement each. Note that, as 

Bending and Torsion customary, displacements are measured relative 

to the elastic axis (Fig. 2). Each clamped beam 
is assumed to be discretized by the Galerkin method in conjunction with two shape functions per 
displacement component. For bending, the shape functions are chosen as the eigenfunctions of a 
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uniform cantilever beam 


Kir = sin ft r Xi - sinh (5 r Xi - sin ^ r ^ + Sin ^ (cos (3 r Xi - cosh f3 r Xi), r = 1,2; i 

cos jj r Li + cosh o r Li 


f,w,e 

(105) 


and for torsion, the eigenfunctions of a uniform clampe-free shaft 

<j>i>ir = sin(2r - l)nXi/2Li, r — 1,2; i- f,w,e (106) 

where Li is the length of the cantilever beam. In this regard, we note from Fig. 2 that the fuselage 
is modeled as two cantilever beams clamped at Of, one pointing to the aircraft nose and the other 
to the tail. New quantities entering into the first-order equation are Cj- 1 * and E ( f l) , which can 
be obtained from C / and Ej by using 0 f = 0 { p + 9 { j ] and letting the components of Of be 
small, F (1) , M (1) , Q ui and (i = f,w,e), which are given by Eqs. (54), and $ ui and 
which contain <t> uir and Kir as given above. Moreover, the stiffness matrices are obtained from the 
potential energy as follows: 

V = \ f Df [ Elfz (d 2 Uf y /dx )) 2 + EIfy (d 2 u fz /dx 2 f ) 2 + GJ fx &ip fx /dx f ) 2 ] dD f 

+ f Dw \EI W (3Pu wz /dx 2 w ) 2 + GJ W {d^ wx /dx w f\ cLD w + J De [EI e {cEu ez /dx 2 e f 
+GJ e [dEex/dXef] dD e = \Y,x(<£,.i K ux<\ux + qJiA^q^) (107) 

where 

K., = f D (Kf) T di*s[EIn EI, y )K,dD„ K ui = / D £/ f (4>" 

( 108 ) 

K i>i = J D GJi(^) T %idD u i = f,w,e 

are the desired stiffness matrices, in which primes denote differentiations with respect to Xi. 

Equations (63) for the model in question are of order 76, but the equations are linear. However, 
in the case of certain aircraft maneuvers, the systems is time-varying. 

To demonstrate the ideas, we consider two cases, steady level flight and steady level turn ma- 
neuver. 


i. Steady level flight 

For steady level flight, the zero-order velocities are defined by 

v (°) = C (0)[y(0) 0 0 jT = constant) w <°) = 0 (109) 

where is the aircraft forward velocity. From Eqs. (49), we conclude that the zero-order mo- 
menta are 

Py| = mV ( p = constant, p^J = S <0) V^ 0) = constant (110) 

Hence, from the second line of Eqs. (43), we have 

F (0) = 0, M (0) = 0 (111) 
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The implication of Eqs. (Ill) is that, for steady level flight, the forces and moments due to the 
engine thrust, aerodynamic forces, gravitational forces and control forces balance out to zero. The 
angle of attack can be expressed as 

ctf = tan ~\vf*/V$) = a Q = constant (112) 

For level flight, we have V ;(0) = ( P {0) — 0, so that the pitch angle is equal to the angle of attack, or 

0<°> = af (113) 

Moreover, because vj^ = 0, the sideslip angle is zero, 

(3f 0) = tan -'(Vff/Vjg) = 0 (114) 

In view of this, and due to the symmetry of the gravitational and aerodynamic forces, the side 
force Fy 0} and the roll and yaw moments, Mx 0> and are automatically zero. We assume 
that V (0) = 5000 in/s and consider a flight at a 25000 ft altitude, so that the speed of sound is 
1016.1 ft/s and, hence, the Mach number is 5000/(1016.1 x 12) = 0.41. From Eq. (109), we have 
Vj?) = 5OOO[cos0( o) 0 sin0 (o) ] 7 \ Then, using Eqs. (44) in conjunction with Eqs. (45)-(48), Eqs. 
(Ill) yield 

Fj 0) = 2 F^ 0) - 695.3237 cos 2 0 (o) - 12973 sin 0 {o) + 2149.5358^ 0) cos 0 (o) sin 0 (o) 

+ sin 2 0 (O) ( 164224.6592 + 2149.5358^ 0) tan0 (o) + 164919.9829 tan 2 0 (o) ) = 0 
Fj 0) = - 2149.5358^ 0) cos 2 0 (o) + cos 0 (o) { 12973 - 202450.4921 sin0 (o) ) 

+ sin 2 0 (o) (-2149.5358^ o) - 239285.6777 tan 0 (o) - 36835.1856 tan 3 0 (o) ) = 0 

(115) 

M( 0) = - 27.92 F^ 0) + (-16069.8929 - 553904.5798tf< 0) ) cos 2 0 (o) 

-52019.1354 sin 0 (o) + (-2.0058 x 10 6 - 116986.5503^ 0) ) cos 9 {0) sin 6 {0) 

+ sin 2 0 (o) (3.0203 x 10 6 - 553904.5798^ o) + (-1.6634 x 10 6 
- 116986.5503^ o) ) tan0 (o) + 3.0364 x 10 6 tan 2 0 (o) + 342353.8625 tan 3 0 (o) ) = 0 

which can be solved for the pitch angle 9 i0 \ the engine thrust F^’ and the elevator angle <5 f <0i . Solv- 
ing the nonlinear equations (115), we obtain = 0.0667 rad, F^ = 431.6465 lb and = 

—0.2703 rad, so that the zero-order control vector is given by = [F^ 0) 0 Se 0> 0] r = 

[431.6465 0 — 0.2703 0] r . Hence, the control force vector can be written in the matrix form 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

2 0 143.4894 0 

0 0 0 -1905.7803 

0 0 -2149.5358 0 

0 803506.9168 0 -153855.5426 

-27.92 0 -561713.8624 0 

0 53637.0429 0 529190.6764 

= [0 0 0 0 0 0 824.5038 0 581.0801 0 139795.5057 0] T 


p(0) _ £(0) u (0) _ 
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In the case of steady level flight, the aircraft experiences static deformations due to zero-order 
forces. Because static deformations are constant and deformations are first-order quantities, we 
denote all quantities involved by a subscript c and a superscript (1). Consistent with the zero-order 
results and using the first of Eqs. (50), we write 

vf c = cf {[V (1) 0 0] T - u >f e = 0 ; s uic = 0 , s^ c = 0 ,i = f,w,e (117) 

where VO) j s the first-order forward velocity. Then, from Eqs. (59), the corresponding momenta 
are 


p"/c = mv£>, = S««Vf + S< 0 ’V<" 

and, from Eqs. (50), we conclude that 

F< 1} = 0, M< 1} = 0 

KuiQuic + Q uic — 0, “t - 0, i f i ^ 


(118) 


(119) 


which represent algebraic equations to be solved for the first-order pitch angle, the first-order 
elevator incident angle, the first-order engine thrust and the static generalized displacements q^c 
and q $ ic . We note that = [Fxl } 0 F { z l ] \ T , = [0 M$ 0] r , Q uic and Q VJC are all 

functions of q uic and q^ c (i = f,w,e). 

The stiffness matrices are as follows: 


k f - 

A u/ — 




Kf= 


8687.0713 

-3181.9271 

0 

0 


-3181.9271 

293633.2239 

0 

0 


0 

0 


12883.4627 
-4643.2903 

8891.2109 -28248.0864 0 

-28248.0864 231517.6212 0 


0 

0 

-4643.2903 

419984.1291 

0 

0 


0 

0 

1.5462 

1.0134 


0 

0 


1.0134 

9.0414 


x 10 8 , K$ f = 


13097.4380 -36822.7531 
-36822.7531 363078.6024 

1.0781 1.9339 

1.9339 6.7838 


10 8 


jsR _ jsL _ 
uw uw 


tsR tsL 
n ue ue 


581.3939 

-1488.0598 ' 

K* 

=k l = 

' 2.2974 

3.6438 

-1488.0598 

13050.3375 


t \)W 

3.6438 

13.5871 


( 120 ) 


x 10 7 


406.2228 

-696.8287 ‘ 

, K 1 } 

=k l = 

' 1.5455 

1.4924 

-696.8287 

10685.4900 


“i/»e 

_ 1.4924 

10.0858 


10 6 


Ke = 


2126.9677 -4303.4586 
-4303.4586 52929.3173 


Ke = 


5.0679 

5.3834 


5.3834 

31.9145 


10 6 


where the superscripts F, A, R, L and V denote the fore part, aft part, right half, left half and verti- 
cal (stabilizer), respectively. Assuming that V il) = -5 in/s, so that = V^ 1 ) [cos 6 |(0) 0 sin 0 ( F] T 
+ V’(o)[-0(i)s0(°) o e^cd^f = [-4.6953 - 333.02720^ 0 - 4.7316 + 4988.897O0 (1) ] r , we 
can solve Eqs. (119) and obtain 


d w = - 0.00088 rad, F^ = -3.1169 lb, <$< l > = -0.0047 rad 
q£ /c = [0 0 0.1206 - 0.0015] r in, qj /c = [0 0] T rad 
q£ /c = [0 0 0.0585 0.0038] r in, q$ /c = [0 0] r rad 
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qLc = QLc = [-3.4537 - 0.2918] r in, cg wc = q^ c = [0.0021 - 0.0002] r rad 
qL = quec = [-0.0279 0.0006 ] t in, qj ec = qj ec = [0.0024 - 0.0002] r rad 
qL = [0 0] r in, qj ec = [0 Of rad 

Assuming that the first-order state is measured from the static elastic displacements posi- 
tion, the first-order state equations for steady level flight can be written in the customary form, Eq. 
(63), where A and B are constant coefficient matrices. This requires the mass matrices, defined by 
Eqs. (15). They can be computed as follows: 

Mn = 

0 131.1002 0 1 To 0 -4.2385 ] 

Mio = S T = -131.1002 o 0 + 0 0 0 qF 


0 

131.1002 

0 " 


0 

0 

-4.2385 " 

-131.1002 

0 

0 

+ 

0 

0 

0 

0 

0 

0 


4.2385 

0 

0 


0 

0 

3.2026 " 

r o 

4.2385 

0 1 

0 

0 

0 

Qufy2 + - 4 ' 2385 

0 

Qufzl + 

3.2026 

0 

0 

L 0 

0 

0 


0 

0.0807 

0.0128 


-0.0807 

0 

0 


-0.0128 

0 

0 


r o 

0 

-0.1252 " 

Que2 + 0 

0 

0 

0.1252 

0 

0 


0 0.0219 


<3 = 


[ -0.0219 0 

0 0 
4.2385 -3.2026 


0 

0 

4.2385 


0 

0 

-3.2026 


Moo — J — 


182862.0586 

4.7745 

-37732.0967 “ 


0 

-863.2814 

0 

4.7745 

566132.2704 

81.1583 

+ 

-863.2814 

0 

-21.0393 

-37732.0967 

81.1583 

704093.8389 


0 

-21.0393 

0 


0 

367.2116 

0 

T 42.0786 

0 

-863.2814 " 

367.2116 

0 

7.3021 

Qufy2+ 0 

42.0786 

0 

0 

7.3021 

0 

1V [ -863.2814 

0 

0 


Qufzl + --- < 122 > 


6.9600 

3.2162 

-20.3065 ‘ 


0 

35.9911 

0 

3.2162 

7.5428 

-1.2425 

Que 2 + 

35.9911 

0 

12.8989 

20.3065 

-1.2425 

-0.5828 

0 

12.8989 

0 


0 

-4.0598 

0 


-4.0598 

0 

4.0183 


0 

4.0183 

0 


" -21.0393 

7.3021 

0 

0 ] 

' 0 

0 

2.4862 

-0.7090 1 

0 

0 

-863.2814 

367.2116 + 

0 

0 

0 

+ 

tH 

3 

U- 

o 

863.2814 

-367.2116 

0 

0 

0 

0 

0 

0 J 


0.7090 -2.3655 0 0 


K ; = 


68.4129 

0 


72.7595 

0 


-344.8056 -402.1912 


M 33 = 


2.2363 -1.5318 

-1.5318 8.4942 

0 0 

0 0 


0 

0 

2.0922 

-1.5310 


0 

0 

-1.5310 

7.9257 
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M& = 


17.3669 

10.8317 


10.8317 

32.5338 


K = 


34.1002 

32.0129 


32.0129 

65.4667 


Before we can compute A and B, we must still generate the damping matrices C ui and ( i = 
f,w,e). To this end, we assume that the damping functions Cui and ( i = f,w,e ) are all 

constant, so that the damping matrices, Eqs. (19), reduce to 


B'ui — ^ui^ui 5 ^ V 1 ? — 7 , ) /j 1C, e 


(123) 


Equations (123) state that the damping matrices are proportional to the stiffness matrices, which 
permits us to write the relations 37 

c ui = 2C/A if, = 2C/AJ,; 2 , i = f,w,e (124) 


where £ is a structural damping factor and A^ 2 and A^ 2 (i = /, u\ e) are the lowest natural 
frequencies of the respective components. We assume that £ = 0.03. Moreover, we obtain the 
component natural frequencies by solving the eigenvalue problems 

det \K[, - Af,M£] = 0, det \K* - A 1,M*\ = 0, det[JT* - Aj^M*] = 0 
detyct - aLmZ,] = 0 , .... det [K* - A^A/£j = 0 , det [*£ - Aj,A/£] =0 

with the results 


yfc = 59.1046 rad/s, = 54.1788 rad/s, ^A^ = ^/a£ = 36.2911 rad/s 

= v'AS = 72.1276 rad/s, = 131.9557 rad/s, (126) 


= 221.8479 rad/s, y/X^ = 61.9940 rad/s, y] X = 253.9293 rad/s 
^ = yfc = 294.1938 rad/s y / X^ = 384.8014 rad/s 
Hence, using Eqs. (120), (123), (124) and (126), the damping matrices are 



8.8187 -3.2301 0 0 

-3.2301 298.0817 0 0 

0 0 13.0786 -4.7136 

0 0 -4.7136 426.3467 



9.8465 -31.2832 0 0 

-31.2832 256.3929 0 0 

0 0 14.5047 -40.7792 

0 0 -40.7792 402.0894 


r F = I 41818 - 3661 27408.5814 

W/ - 27408.5814 244530.1083 


r A _ 104342.4068 187172.4172 

W/ 187172.4172 656562.6583 


/~<R s^L 

^ UW ^ U W 


0.9612 -2.4602 
-2.4602 21.5761 


r R r L _ 5428.3762 8609.7859 

8609.7859 32104.4936 


fyR f' iL 

^ue ^ue 


0.3379 -0.5797 c r _ c l 
-0.5797 8.8888 ’ We We 


315.1931 304.3728 

304.3728 2056.9682 



0.9671 -1.9568 

-1.9568 24.0669 



790.2110 839.4058 

839.4058 4976.2592 


( 127 ) 
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Then, the coefficient matrices in Eq. (63) can be shown to be 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

-333.0272 

0 

5000 . . . 

-0.0014 

-0.0001 


0 

0 

0 

0 

-5000 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

36.3883 

19.0770 

0 

-0.0212 

-0.0020 


0 

0 

0 

40.8390 

41.4025 

0 

0.1414 

0.0257 


0 

0 

0 

455.5429 

0 

0 

-31.2223 

-9.6060 


0 

0 

0 

538.2891 

0 

0 

63.7274 

-136.8998 


0 

0 


0 

0 




(128) 

0 

0 


0 

0 




0 

0 


0 

0 




2 

0 

143.4894 

0 




0 

0 


0 -1905.7803 




0 

0 


0 -3097.7273 





Moreover, the feedback control vector is 

„<■> = [F™ <5™ <S<‘> ^ 1 >] T (129) 


in which it is assumed that the right and the left ailerons rotate by angles of the same magnitude <5£ x) 
but of opposite sense and that the right and left elevators rotate by angles of the same magnitude 
and sense. 


Choosing the weighting matrices Q and R in the performance index, Eq. (66), as follows 


1 

0 

0 

0 

0 

0 

0 . 

. 0 ' 


0 

1 

0 

0 

0 

0 

0 . 

. 0 


0 

0 

1 

0 

0 

0 

0 . 

. 0 

r 

0 

0 

0 

1 

0 

0 

0 . 

. 0 


0 

0 

0 

0 

1 

0 

0 . 

. 0 

R = 

0 

0 

0 

0 

0 

1 

0 . 

. 0 


0 

0 

0 

0 

0 

0 

0 . 

. 0 


0 

0 

0 

0 

0 

0 

0 .. 

. 0 



10 

0 

0 

0 


0 

5 x 10 8 
0 
0 


0 

0 

10 8 

0 


0 

0 

0 

10 8 


(130) 


solving the steady-state Riccati equation, Eq. (73), and using Eq. (68), we obtain the gain matrix 


0.2994 

0 

0 

0 

0 

0 

0 

0 

-0.1016 

0 

0.0001 

0 

0.0003 

0.5222 

0 

0.1024 

-1112.0150 

0 

-1.8943 

0 

-0.0187 

1.9350 

0 

0.5103 

-0.0006 

0 

0 

0 

0.0001 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


( 131 ) 
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Then, solving the closed-loop eigenvalue problem, Eq. (75), we obtain the closed-loop eigenvalues 

A u = -0.1038 ± 0.0868 i, A 3 , 4 = -0.1174 ± 0.2033 i, A 5 = -0.2349, 

A 6(7 = -0.2918 ± 0.3064 i, ... , A 69i70 - -40.5617 ± 599.4532 i, 

A 71i72 = -40.7342 ± 604.4685 i, A 73>74 = -79.0030 ± 791.6011*, (132) 

A 75) 76 = -211.4842 ± 1063.7436 i 


Clearly, all the eigenvalues are real and negative or complex with negative real part, so that the 
closed-loop first-order system is asymptotically stable. The implication is that any disturbances 
from the steady level flight are driven to zero. This is in contrast with the open-loop eigenvalues, 
the eigenvalues of A, the first four of which are zero and the fifth is real and positive. 

Finally we consider the response of a closed-loop system to a gust acting on the wing and 
having the linearly distributed form 

t) = [ 0 0 - 0.5(3 + x*/L w ) sin7rt~(l - t)] T , 0 < x* < L w 
fwi x wi t) = [° 0 - 0.5(3 - x^/L w ) sin TtU(\ - t)] T , 0 < < L w 

where sin 7r^(l — t) represents a half-sine pulse, in which *(1 — t) is a rectangular function 
of unit amplitude and unit length. Inserting Eq. (133) into Eqs. (54), we obtain the generalized 
force components of the disturbance vector F ex t entering into Eq. (74), which can be integrated 
to obtain the system response. Figures 3-5 show the response for the rigid body variables and a 
selected number of elastic variables, and Fig. 6 shows the control inputs. 





Figure 3. Rigid Body Displacements 
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Figure 6. Control Inputs 

Next, we turn our attention to the observer. To this end, we assume that and 0 { j } are 
available as part of the output of an inertial navigation system and that there are four sensors 
measuring velocities at the points P*. (k = 1, 2, 3, 4) of each aircraft component. The coordinates 
Xk, Uk of Pfc in component body axes, together with the type of velocities measured, are listed in 
Table 1 . Moreover, the sensors locations are shown in Fig. A. 


Table 1 



Fore 

Fuselage 

Aft 

Fuselage 

Right 

Wing 

Left 

Wing 

Right 

Horizontal 

Stabilizer 

Left 

Horizontal 

Stabilizer 

Vertical 

Stabilizer 

Pi -21,2/1 
Velocities 

111.21,36.25 

111.92, 35.0 
V^andV), 1 ’ 

120.35,58.23 

y-t 1 ) 

V WZ 

120.35, -58.23 
y(l) 

V WZ 

50.54, 32.0 

yd) 

V €Z 

50.54, -32.0 

50.0, 42.0 

v£> 

P2 • 22,1/2 
Velocities 

111.21, -36.25 

yW 

111.92, -35.0 

W 

120.35, -28.0 

y(!) 

v wz 

120.35, 28.0 

yO) 

v wz 

46.6, -20.0 

yd) 

v ez 

46.6, 20.0 

yd) 

v ez 

41.0, -28.0 

P3 ■ 23 , J/3 

Velocities 

277.21,24.0 

279.79, 25.0 
Vf? and ^’ 

327.66, 25.07 

y(!) 

v wz 

327.66, -25.07 

y(D 

v wz 

127.0, 19.0 

v£> 

127.0, -19.0 

yd) 

v ez 

112.0.27.0 

W 

P 4 : 24 , 2/4 
Velocities 

277.21, -24.0 

279.79, -25.0 

327.66,-10.0 

y(i) 

v wz 

327.66. 10.0 

yd) 

V WZ 

127.0, -12.0 

127.0, 12.0 

y(D 

1 151 

113.0, -16.0 
v * 11 

v ez 


To compute the matrix C, Eq. (103), relating the output vector to the state vector, we use Eq. (100) 
and write 


0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 

1411.5402 

- 1066.5377 

0 

0 ... 

0 

0 

0 

21145.5039 

0 ... 

0 

- 21145.5039 

15977.2118 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 
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( 134 ) 


0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 

4.2544 

41.7049 

- 7.3013 

0 ... 

0 

. . . — 402.3892 

0 

0 

0 ... 

0 

. . . - 63.7322 

- 624.7574 

109.3773 

0 ... 

0 

0 

0 

0 

0 ... 

0 

0 

0 

0 

0 ... 

0 


Then, inserting Eq. (134) into Eq. (103), we obtain 





" 1 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

1 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0.6890 

- 0.5206 

0.0026 

0.0019 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 1.8548 

1.8590 

1.4015 

- 1.4046 

4.1411 

4.1328 

- 3.1290 

- 3.1227 


0 

0 

0 

0 

0 

0 

1.5244 

- 1.1518 

- 0.0066 

0.0050 


0 

0 

0 

0 

0 

0 

- 2.4135 

1.8236 

- 0,1611 

0.1217 


_ 0 

0 

0 

0 

0 

0 

1.8189 

- 1.3970 

0.1101 

- 0.0832 


0 

0 

0 

0 

0 

0 

0 

0 

- 0.00003 

0.0123 . . . 

- 0.0002 

0 

0.0008 

0.0320 . . . 

- 0.0004 

0 

0.0011 

- 0.0373 . . . 

0.0004 

0 

- 0.0015 

0.0871 . . . 

- 0.4412 

- 0.6612 

0.0002 

- 0.0260 . . . 

2.1761 

- 1.4006 

0.0009 

- 0.0385 . . . 

- 1.4522 

0.9932 


( 135 ) 


We assume that the excitation noise v(£) and observation noise w(t) represent white noise pro- 
cesses with intensities V and IV, respectively, and choose 


V = diag[ 100000 100000 100000 100000 100000 100000 10~ 9 10~ 9 

10-9 10 -9 10 -9 10 -9 10 -9 10 -9 10 -9 jq -9 jq -9 10~ 9 10~ 9 10~ 9 10~ 9 10~ 9 

l 0 -9 1q -9 1q -9 1q -9 1q -9 1q -9 1q -9 jq -9 jq -9 10 - 9 1(r 9 10~ 9 10~ 9 10~ 9 0 ... 0 ] 

(136) 


W = diag[10~ 12 10~ 12 10~ 12 10~ 6 10~ 6 10~ 6 1 ... l] (137) 

so that, from Eq. (91), the steady-state Riccati equation 

AQ + QA t + V - QC t W~ x CQ = 0 (138) 


yields 

Q = diag [ 0.0003 0.0003 0.0003 0.3162 0.3162 0.3162 ... 0 0 ] 


(139) 
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Hence, from Eq. (90), the optimal observer gain matrix is 

K 0 = QC T W~ l 


1623 x 10 8 

0 

0 

0 

0 

0 

0 

0 

3.1623 x 10 8 

0 

- 0.0003 

0 

0.0050 

0 

0 

0 

3.1623 x 10 8 

0 

- 0.0050 

0 

0 

0 

- 332.6945 

0 

316227.7660 

0 

0 

0 

0 

0 

- 4995.0051 

0 

316227.7661 

0 

0 

0 

4995.0050 

0 

0 

0 

316227.7661 . . . 

0 

0.0018 

- 0.6163 

0.0023 

- 0.0471 

- 0.0003 

0.1358 

0 

- 0.0013 

0.4656 

- 0.0018 

0.0348 

0.0002 

- 0.1029 

0 

3.4719 

0.0042 

4.5563 

0.0001 

- 0.5724 

0.0003 

0 

- 2.6233 

- 0.0032 

- 3.4428 

- 0.0001 

0.4346 

- 0.0002 

0 

0.0016 

- 0.5422 

0.0021 

- 0.0414 

- 0.0003 

0.1193 

0 

0 

- 0.0014 

0 

0 

0 

0 

0 

0 

- 0.0001 

0 

0 

0 

0 

0 


(140) 


Finally, solving the eigenvalue problem for A — K 0 C, we obtain the observer eigenvalues 

Ai = -0.0103, A 2 , 3 = -0.7807 ± 1.6618z, A 4)5 = -0.1874 ± 2.2166?, 

A 6 = -5.7278, A 7 , 8 = -5.7237 ± 36.9798?, A 9 , 10 = -3.7696 ± 45.6080?, 

An, is = -4.6367 ± 62.6348?, Ai 3 ,i 4 = -4.2692 ± 68.7929?, . . . , 

A 63 ,64 = -40.5621 ± 599.4494*, A 65 , 66 = -40.7343 ± 604.4687?', 041 } 

A 6 7,68 = -79.0066 ± 791.5821?, A 69)70 = -211.4807 ± 1063.7432?, 

A 71 . 72.73 = —316228, A 7 4 , 7 5 , 7 g = —3.1623 x 10 8 

The performance of the observer design can be demonstrated by simulating the response of the 
combined system, defined by Eq. (84), to an initial state and initial observer error. To this end, we 
choose the values 

x (1) (0) = [5 5 5 0.005 0.005 0.005 0.2 0.1 0.2 0.1 0.2 0.1 -0.2 -0.1 -0.2 

-0.1 -0.2 -0.1 -0.2 -0.1 -0.2 -0.1 -0.2 -0.1 0.05 0.01 -0.05 
-0.01 0.05 0.01 -0.05 -0.01 0.05 0.01 -0.05 -0.01 0.05 0.01 0... 0] r 
e(0)= -0.15x (1) (0) 

Figures 7-9 show the response for a selected number of rigid body and elastic variables 
observer estimates. Figure 10 shows the control inputs as given by Eq. (82). 


(142) 
and their 
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Estimated 


Figure 9. Generalized Elastic Displacements 




ii. Level steady turn maneuver 

We consider the case in which in the zero-order problem the aircraft flies at a constant velocity 
around a circular path of radius R in the horizontal Y-plane. In this case, it is convenient to 
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refer the rigid body motions to a set of axes x\y\Z\ obtained through a rotation p) (0) about Z, where 

= f l = constant. It is not difficult to see that axes Xi,yi and z\ represent a set of cylindrical 

- ( 0 ) 

axes t, n and Z, where t is tangent to the circle and n is normal to it. Denoting by Ry the velocity 
of Of in terms of cylindrical components, the kinematical relation corresponding to the first of 
Eqs. (43) can be written as 


where 


A(0) = [< 


R 


R 


(0) 

fn 


Z {0) } T = [RCl 0 Of = Cf )T 


y(0) 



c0<°> 0 -s0 (o) 

s9^s<p^ ccj)^ c9^s<j>^ 

s9^c 4>^ — s <f>^ c9^c<j>^ 


(143) 


(144) 


is the matrix of direction cosines between tnZ and the fuselage body axes XfyfZf, obtained from 
Cf, the first of Eqs. (2), by letting tp = 0 and replacing 6 and 4> by 9 {0) and f 0) , respectively. 
Similarly, the second of Eqs. (43) can be written as 

ef = [ 0 (o) 9 {0) ^ (0) ] = [0 0 fl] r = {E ( p)- l uf (145) 

where can be obtained from Ej, the second of Eqs. (2), by replacing 9 and <j> by 9 {Q) and 4> {0 \ 
respectively. Equations (143) and (145) yield 

V { p = [vP Vp VpY = C ( P Ry 0> = RCl{c9 (0) s0 (O) s0 (O) s9 (0) c<j> (0) } T = constant (146) 
and 

u>P = [u>P u>P ^p] T = Ep6p = fl[— s<9 (0 ^ c# (0 ^s<A (0) c9^s<f>^} T = constant (147) 
so that, from Eqs. (49), we have 

Pyj = mVp + S (0>T u) ( p = constant, ppj = S (0 ^Vp + J^ujp = constant (148) 
It follows that the equations of motion, the last two of Eqs. (43), reduce to 


-Cbfp { v) + F (0) = 0, -V} 0 VP f - u>J 0) pi°j 4- M<°> = O (149) 

which are independent of time. 

To determine the parameters defining the steady level turn maneuver, we choose the turn radius 
R and angular velocity Cl and solve Eqs. (149) for the bank angle <^°\ pitch angle 9 (Q) and control 
vector = [F e (0) Sa 0> Se 0) Sl 0> ] T . Hence, assuming the values R = 1.5 mi = 95037 in and 
Cl = 0.0526 rad/s, so that RCl = 5000 in/s, and using Eqs. (146) and (147), we have 



’ C0<°> ' 


* -S0<°> ‘ 

VP = 5000 

sO^scf)^ 

, jp = 0.0526 

c9^s<j>^ 



c0<°)c</><°) 


We consider a flight at a 25000 ft altitude, so that the speed of sound is 1016.1 ft/s and, hence, the 
Mach number for the flight is 5000/(1016.1 x 12) = 0.41. Inserting Eqs. (150) into Eqs. (149) in 
conjunction with Eqs. (148) and solving the resulting transcendental equations, we obtain 

6> (0) = 0.0986 rad, <£ (0) = 0.6160 rad 

(151) 

Fjp = 468.74291b, <^ 0) = -0.0028 rad, S ( e 0) = -0.3233 rad, ^ 0) = -0.3445 rad 
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Then, the zero-order control force vector can be written as 


p(0) = g(0) n (0) = 


0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

- 27.92 

0 


0 

0 

0 

0 

0 

0 

1.0370 

0 

- 13.7394 
1.0746 x 10 6 
- 164.5690 
65027.5713 


0 

0 

0 

0 

0 

0 

178.2611 

0 

- 2875.1495 

I . 4415 

- 750586.5476 

II . 7592 


0 

0 

0 

0 

0 

0 

77.3023 

- 1916.1116 

0 

- 154689.6019 

- 6240.6944 

532059.4480 


468.7429 

- 0.0028 

- 0.3233 

- 0.3445 


0 

0 

0 

0 

0 

0 

844.5606 

654.9960 

692.9931 

50641.1455 

171482.9784 

- 182061.4664 


( 152 ) 


As in the case of steady level flight, the aircraft experiences static deformations in the steady 
level turn maneuver as well. Denoting the corresponding constant quantities by the superscript 
(1) and subscript c, using the first of Eqs. (50) and letting V (1) be the cylindrical coordinates 

-U) 

counterpart of K f , we can write 


= C} 0) [V (1) - Cf T Vf] = constant, ujf c = 0 
0, i /*, tc, c 


(153) 


The momenta are all constant and can be expressed in terms of the zero-order and first-order 
velocities and the static deformations using Eqs. (59). 

A constant solution of the first-order equations, Eqs. (50), can be obtained by letting the left 
sides be equal to zero and solving for 4>c\ 0^, uf = [FjV 6$ f and the static 

deformations q uic and q^ c (i = f, w, e). Assuming that V (1 ) = [-5 0 Of in/s and using the first 
of Eqs. (153), we obtain 

-4.9754 - 492.12170c 1 ) 

-0.2843 + 2874.81210c 1} + 486.46524 1} (154) 

-0.4017 + 4061.19090c 1 ) - 344.35624 1 ) _ 

so that Eqs. (50) yield 

4 !) - -0.0011 rad, 0f = -0.0013 rad 

uf = [-4.0727 - 0.0002 - 0.0063 O.OOOlf 

q£ /c = [-0.0005 0.0001 0.1458 - 0.0019f , qj /c = [-0.0001 Of 

q£ /c = [-0.0408 - 0.0008 0.0827 0.0039f , q} fc = [0.0003 - 0.0001] T 

<£* = [-4.1523 -0.351lf,qJ BC = [0.0025 - 0.0003f 

qLo = [-4.1765 - 0.3526f, q^ c = [-0.0025 0.0003f 
q« ec = [-0.0473 0.0001] T , q* c = [0.0029 - 0.0002] r 
q«ec = [-0.0489 0.00008f , cfi ec = [-0.0029 0.0002f 
q^ec = [0.0032 0.0014] r , q£ ec = [0.0017 - 0.0001] T 
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The mass matrices, Eqs. (15), for the steady turn case are as follows: 


M 12 - = 


0 

130.4114 

0.1477 " 


0 

0 

- 4.2385 " 

- 130.4114 

0 

0 

+ 

0 

0 

0 

- 0.1477 

0 

0 


4.2385 

0 

0 



0 

0 

3.2026 " 


0 

4.2385 

0 " 

+ 

0 

0 

0 

Qufy2 “1“ 

0 

0 

0 


- 3.2026 

0 

0 

- 4.2385 

0 

0 


F 

Qufyl 


Qufzl + 


+ 


0 

- 0.1095 

- 0.0173 


0.1095 

0 

0 


0.0173 

0 

0 


L i 

Que 2 + 


0 

0 

0.1252 


- 0.1252 

0 

0 


V 

Quel 



0 

0 

0.0219 " 


0 

0 

0 

0 

+ 

0 

0 

0 

& 2 . <3 = 

4.2385 

- 3.2026 

0 

0 


- 0.0219 

0 

0 

0 

0 

4.2385 

- 3.2026 


hi 22 — J — 


182795.8456 

- 18.6695 

- 37747.5473 " 


- 0.0025 

- 863.2814 

0 

- 18.6695 

566091.8052 

74.4042 

+ 

- 863.2814 

0 

- 21.1022 

- 37747.5473 

74.4042 

704068.0911 


0 

- 21.1022 

- 0.0025 


F 

Qufyl 


+ 


+ 


+ 


0.0010 

367.2116 

0 


42.2044 

0 

- 863.2814 " 

367.2116 

0 

7.3208 

Qufy2 "I” 

0 

42.2044 

0.0012 

0 

7.3208 

0.0010 

- 863.2814 

0.0012 

0 


Qufzl + 


6.9557 

3.2162 

- 20.3065 " 


’ - 0.0075 

35.9911 

0 

3.2162 

7.5394 

- 1.2456 

Que 2 + 

35.9911 

0 

12.8965 

- 20.3065 

- 1.2456 

- 0.5837 

0 

12.8965 

- 0.0075 


( 156 ) 


v 

Quel 


0.0022 - 4.0598 0 

- 4.0598 0 4.0187 

0 4.0187 0.0022 
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" - 21.1022 

7.3208 

- 0.0012 

0.0005 


’ 0 

0 

2.4362 

- 0.7090 " 

— 

0 

0 

- 863.2814 

367.2116 

+ 

0 

0 

0 

0 


863.2814 

- 367.2116 

0 

0 


0 

0 

0 

0 


Qufyl 
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+ 


2.4862 

- 0.7090 

0 

0 

17.3669 

10.8317 


" 0.7090 

- 2.3655 

0 

0 " 


ii 

68.3900 

72.7325 

0 

0 

0 

0 

Qufz2'> * * 

0 

0 

0 

0 

0 

0 

- 344.8056 

- 402.1912 


- 0.7090 

2.3655 

0 

0 

10.8317 

32.5338 


0 

0 

2.4862 

- 0.7090 


, K l 


0 

0 

- 0.7090 

2.3655 

34.1002 

32.0129 


32.0129 

65.4667 


Next, we consider the time-varying part of the first-order problem, Eqs. (50) with all quantities 
measured from the constant static solution. To this end, we use Eq. (63) in which the coefficient 


49 



matrices are given by 


0 

0 

0 

0 

0 

0 

0.0004 

- 0.0001 

0 

0 

0 

0 

- 596.0114 

0 5000 

. . . - 0.0008 

- 0.0012 

- 0.0001 

0 

0 

0 

0 

-5000 

0 

0.0010 

- 0.0008 

0 

0 

0 

0 

0 

0.0529 

0 

0 

0 

0 

0 

0 

0 

- 0.0524 

0 

0 

0 

0 

0 

0 

0 

0 

174.9021 

42.3914 

0 

... - 76.1004 0.1415 

0.0257 

0 

0 

0 

370.8328 

- 25.9627 

0 

0.0604 

- 31.2428 

- 9.6184 

0 

0 

0 

438.1920 

- 30.6787 

0 

0.1026 

63.7171 

- 136.9109 


0 


0 

0 


0 




0 


0 

0 


0 




2 


1.0326 

176.1982 


104.3594 




0 


0 

0 


- 1901.4388 




0 


- 10.2357 

- 2143.1016 

0 




0 


800910.3941 

1.8988 


- 153505.0524 



— 

27.92 

- 116.9527 

- 561836.0209 

- 8425.0411 




0 


64654.3829 

11.7819 


527985.1550 




0 


0 

0 


- 7461.0801 




0 


0 

0 


- 3090.6705 





and in which the feedback control vector has the form u (1) = [F { ^ ] 8^ 5 ( e ^ <5r 1, ] T . In the case 
in which u (1) = 0 and F ext = 0, the state equations admit an exponential solution yielding an 
eigenvalue problem. Solving the eigenvalue problem, we conclude that the system is unstable, 
with four eigenvalues being equal to zero and one being real and positive. Using a linear quadratic 
regulator in conjunction with the weighting matrices 


Q = 


10 0 0 
0 10 0 
0 0 10 
0 0 0 1 
0 0 0 0 
0 0 0 0 


0 0 0 0 


0 0 
0 0 
0 0 
0 0 
1 0 
0 1 


0 

0 

0 

0 

0 

0 


0 0 


0 


, R = 


10 0 0 0 

0 6 x 10 8 0 0 

0 0 10 8 0 

0 0 0 10 8 


(158) 


solving the corresponding steady-state Riccati equation, Eq. (73), and using Eq. (68), we obtain 
the gain matrix 


0.2974 

0 

- 0.0163 

0 

- 0.1062 

0 

12.4530 

0.5355 

- 1022.6944 

- 1.0086 

- 71.9234 

1.3797 

- 0.0006 

0 

0.0001 

0 

0 

0 

0 

0 


0 

0 

- 0.0001 

0 

0.0001 

0 

0.0531 

0.1015 

- 1.2660 

- 0.2651 

- 1.4307 

0.3965 

0 

0 

0 

0 

0 

0 

0 

0 


( 159 ) 
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Then, solving the closed-loop eigenvalue problem, Eq. (75), we obtain the closed-loop eigenvalues 

A 1j2 = -0.1031 ± 0.0870*, A 3 = -0.2342, A 4(5 = -0.1194 ± 0.2172*, 

A 6 7 = -0.2998 ± 0.3089*, A 8 9 = -0.7842 ± 1.6226*, 

’ (160) 
A 69> 7o = -40.5623 ± 599.4558*, A 71 , 72 = -40.7352 ± 604.4713* v 

A 73)74 = -79.0053 ± 791.6059*, A 75 , 76 = -211.4850 ± 1063.7465* 

Clearly, all the closed-loop eigenvalues are either real and negative or complex with negative real 
part, so that the closed-loop first-order system is asymptotically stable. Hence, any disturbances 
from the steady level turn maneuver will be driven to zero. 

Finally, we compute the response of the closed-loop system to the gust given by Eqs. (133). 
Figures 11-13 show a selected number of rigid body and elastic variables, and Fig. 14 shows the 
control inputs. 
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Figure 14. Control Inputs 


Next, we turn our attention to the observer. To this end, we must compute the matrix C, Eq. 
( 103), relating the output vector to the state vector. First, we use Eq. (100) and write 


r o ... 

0 

- 0.1811 

0.1369 

0.1282 



o ... 

0 

0 

0 

0.0219 



0 ... 

0 

- 0.0219 

0.0166 

0 



0 ... 

0 

1676.3711 

1275.2629 

- 1242.2537 



0 ... 

0 

3.5685 

- 1.5887 

21090.9413 



0 ... 

0 

- 21090.3030 15935.2420 

0 



0 ... 

0 

0 

0 

0 



( 161 ) 



- 0.0030 

- 0.0054 

0.0009 

0 ... 

0 




- 0.0004 

0 

0 

0 ... 

0 




- 0.0001 

- 0.0006 

0.0001 

0 ... 

0 




27.2581 

51.3897 

- 8.9291 

0 ... 

0 




- 401.1684 

0.3647 

0,1928 

0 ... 

0 




- 63.5214 

- 622.7178 109.2101 

0 ... 

0 




0 

0 

0 

0 ... 

0 



Then, inserting Eq. (161) into Eq. (103), we obtain 


1 

0 

0 

0 

0 

0 

0 

0 
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0 .. 

. 0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 
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- 0.4820 
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- 0.1651 

0 .. 

. 0 

0 

0 

0 

0 

0 

0 

0 

- 2.2285 

1.6958 

5.8787 

- 4.4245 

0 .. 

. 0 

0.0010 

0 

0 

0 

0 

0 

0 

2.2323 

- 1.6989 

2.3810 

- 1.8161 

0 .. 

. 0 

0.0011 

0 

0 

0 

0 

0 

0 

1.4370 

- 1.0831 

0.3731 

- 0.2782 

0 .. 

. 0 

- 0.0015 

0 

0 

0 

0 

0 

0 

- 2.3298 

1.7578 

- 0.5259 

0.3938 

0 .. 

. 0 

0.0002 

0 

0 

0 

0 

0 

0 

1.8150 

- 1.3705 

0.2414 

- 0.1811 

0 .. 

. 0 

0.0009 
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0 


0 


0 


0 

0 

0 

0.0122 . . 

-0.0002 

0 

0.0374 . . 

-0.0004 

0 

-0.0372 . . 

0.0004 

0 

0.0871 . . 

. -0.4412 

-0.6612 

-0.0259 . . 

2.1761 

-1.4006 

-0.0384 . . 

. -1.4522 

0.9932 


(162) 


We assume that the excitation noise v(() and observation noise w (t) represent white noise pro- 
cesses with intensities V and W, respectively, and choose 


V = diag[ 100000 100000 100000 100000 100000 100000 IO -9 1(T 9 10“ 9 IO -9 IO -9 10“ 9 IO -9 

io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 io -9 (163) 

IO" 9 IO" 9 IO -9 10 -9 10 -9 10 -9 IO -9 IO -9 IO -9 0 ... 0] 


and 


W = diag[10 -12 IO -12 IO -12 IO -6 IO -6 IO -6 1 . . . 1] 
so that, from Eq. (91), the steady-state Riccati equation 

AQ + QA t + V - QC T W~ l CQ = 0 

yields 


0 

0.00032 

KJ 

0 

u 

0 

w 

0 

u ... 

0 

w 

0 

... 

0 

V 

0 

0 

0 

0.00032 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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-0.0001 . . . 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0705 

0.0288 . . . 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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(164) 

(165) 


(166) 
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Hence, from Eq. (90), the optimal observer gain matrix is 

K 0 = QC T W~ l 


' 3.1623 x 10 ® 

0 

0 

0 

0 

0 

0 

0 

3.1623 x 10 ® 

0 

- 0.0006 

0 

0.0050 

0 

0 

0 

3.1623 x 10 ® 

0 

0.0050 

0 

0 

0 

- 595.4161 

0 

316227.7660 

0.0003 

0 

0 

0 

0.0001 

- 4995.0050 

0.0003 

316227.7661 

0.0026 

0 

0 

4995.0050 

0 

0 

0.0026 

316227.7661 . . . 

0 

- 0,0427 

- 0.7057 

- 0.4951 

- 0.0554 

- 0.0783 

0.1107 

0 

= 0.0326 

0.5398 

0.3786 

0.0414 

0.0594 

- 0.0839 

0 

3.5739 

- 1.8039 

4.2355 

0.0150 

- 0.4665 

- 0.3316 

0 

- 2.6997 

1.3723 

- 3.1937 

- 0.0107 

0.3542 

0.2518 

0 

0.0003 

- 0.0009 

0.0010 

0 

0 

0 

0 

0.0004 

- 0.0008 

0.0010 

0 

0 

0 

0 

- 0.0001 

- 0.0012 

- 0.0008 

0 

0 

0 

0 

0 

- 0.0001 

0 

0 

0 

0 
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(167) 

Finally, solving the eigenvalue problem for A — K 0 C, we obtain the observer eigenvalues 

Ai = -0.0111, A 2i3 = -0.7858 ± 1.6291*, A 4 , 5 = -0.1935 ± 2.223k, 

A 6 = -5.7234, A 7i8 = -5.7294 ± 36.9785?, A 9 , 10 = -3.7702 ± 45.6096?, 

A u , 12 = -4.6312 ± 62.6280?, . . . , A 63 , 6 4 = -40.5646 ± 599.4514?, 

A 65 ,66 = -40.7353 ± 604.4719?, A 67 , 6 8 = -79.0086 ± 791.5892?, (I68) 

A 69 ,7o = -211.4819 ± 1063.7463?, A 7 i, 72 = -316227.8058 ± 0.0346?, 

A 7 3 = — 316227.8459, A 7 4 t7 5 i7 6 = — 3.1623 x 10^ 

To check the performance of the observer just designed, we simulate the response of the com- 
bined system, Eq. (84) to the initial conditions 

x (1) (0) = [5 5 5 0.02 0.005 0.001 0.2 0.1 0.2 0.1 0.2 0.1 -0.2 -0.1 
-0.2 -0.1 -0.2 -0.1 -0.2 -0.1 -0.2 -0.1 -0.2 
-0.1 0.05 0.01 -0.05 -0.01 0.05 0.01 -0.05 -0.01 (169) 

0.05 0.01 - 0.05 - 0.01 0.05 0.01 0 ... 0] 
e(0) = 0.15x^(0) 

Figures 15-17 show a selected number of rigid body and elastic variables and Fig. 18 shows the 
control inputs as given by Eq. (82). 
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9. Conclusions 


A problem of increasing interest in aeronautics is the simulation of the flight of flexible aircraft 
on a computer. Note that on-line computer simulations can help expedite the design process, and 
possibly reduce the time devoted to flight testing. Moreover, the new class of autonomous aerial ve- 
hicles requires autopilots, which in turn requires real-time simulations. This points to the need for a 
new paradigm in the treatment of flying flexible aircraft, one using the system concept to produce 
a rigorous formulation of the dynamics and control problem and one using potent methodology 
permitting an efficient solution of the problem. This work develops such a paradigm in the form of 
an on-line formulation integrating pertinent material from the disciplines of analytical dynamics, 
structural dynamics, aerodynamics and controls. Moreover, the formulation is cast in a certain 
matrix form ideally suited for computer processing. A perturbation approach permits division of 
the problem into a nonlinear flight dynamics problem for maneuvering quasi-rigid aircraft and a 
linear “extended aeroservoelasticity” problem for the elastic deformations and small perturbations 
in the rigid body translations and rotations, where the solution of the first problem enters as an 
input into the second problem. As a result, there is a different extended aeroservoelasticity prob- 
lem corresponding to each aircraft maneuver. The controls for the flight dynamics problem are 
obtained by an inverse process, which amounts to prescribing a given maneuver and determining 
the controls permitting realization of the given maneuver. On the other hand, the elastic deforma- 
tions and perturbations in the rigid body motions characterizing the extended aeroservoelasticity 
problem are driven to zero by means of feedback controls designed using linear quadratic theory 
in conjunction with a stochastic observer, thus ensuring the stability of the aircraft maneuver. A 
numerical example presents a variety of time simulations of rigid body perturbations and elastic 
deformations about 1 ) a steady level flight and 2) a level steady turn maneuver. 

The integration of the aerodynamics into the unified process is particularly challenging, and 
requires elaboration. For seamless integration, the aerodynamic forces must be referred to the same 
generally noninertial reference frame as that used for all other forces and they must be expressed in 
terms of variables compatible with the variables used throughout the entire formulation. Moreover, 
because the simulation of the system response on a computer is carried out in discrete time, the 
size of the time step is of vital importance. Indeed, on line and real time simulations require that 
the time step be quite small, the time step for the latter being of the order of a minute fraction of 
a second. The implication is that, to be ready to compute the state at the next sampling time, it 
is necessary to be able to compute the aerodynamic forces within the time step; the computation 
of the other forces within the same time step presents no problem. These two requirements are 
quite difficult to satisfy as most aerodynamic theories have been developed for purposes other 
than time response simulations, and the computation of the aerodynamic forces are notorious for 
consuming a great deal of time. In view of this, the development of a new aerodynamic method 
seems highly desirable. Such a method need not be unduly accurate, as robust feedback controls 
should be able to tolerate small errors in the aerodynamic forces. An aerodynamic theory satisfying 
the two requirements outlined above is strip theory. Even though strip theory may not be entirely 
satisfactory for describing the aerodynamic forces acting on whole aircraft, the method is often 
used in aircraft design. In fact it is being used by the same company that provided the data for the 
Numerical Example. The development of an aerodynamic technique suitable for on line, or real 
time response simulations is likely to require a great deal of time and effort. Although the use of a 
more suitable aerodynamic theory is desirable, it is not really necessary at this time. Indeed, at this 
time it is more important to demonstrate how the unified formulation works and how an eventual 
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aerodynamic theory is to be integrated into the overall process. Certainly, such a demonstration can 
provide valuable guidance in the development of an appropriate aerodynamic method to be used 
in conjunction with aircraft time response simulations, thus helping reduce the time and effort 
required for the development in question. 

Finally, it should be pointed out that, with appropriate modifications, the present formulation is 
eminently suited for UAVs, and in particular for autonomous UAV s. The fact that it was applied 
here to an executive jet is due entirely to the ready availability of data from an actual flying aircraft. 
It should also be pointed out that all the time response simulations presented here were carried out 
on a 1 GHz PC using MATHEMATICA. This is particularly important for autonomous UAVs, 
which must be controlled by autopilots, as the required onboard computer is likely to be much 
closer to a PC than to a multiprocessor supercomputer. 
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Appendix - Numerical Values of the Aircraft Parameters 


The values of the pertinent aircraft parameters were provided by an aircraft manufacturer in lumped 
form, and the current formulation assumes distributed parameters. This presents no problem, how- 
ever, as lumped parameters can always be treated as distributed by means of spatial Dirac delta 
functions . 37 

In the first place, we consider the aircraft geometry. To this end, we regard the fuselage as 
consisting of two cantilever beams, a fore part and an aft part, with the origin of both sets of body 
axes at point Of and with axes Xf and xj collinear (Fig. 2), where superscripts denote the fore 
part and aft part. The wing is also divided into two parts, the right half-wing and the left half- 
wing, both with the origin of the respective body axes at Of and Of and with the longitudinal axes 
coinciding with the respective elastic axes. Moreover, the empenage consists of the horizontal 
stabilizer, divided into a right half and a left half, both with the origin of the respective body axes 
at Of and Of, and a vertical stabilizer with the origin of the body axes at Of. The radius vectors 
from 0 / to the corresponding origins are 

>■/» = r/„ = [-5 04 0 38.33 ] t in 

rf e = rf e = [—244.75 0 -43.13] r in, rf e = [-238.97 0 - 24.01] r in 

The formulation calls for matrices of direction cosines between the various component body axes 
and the fuselage body axes, and in particular the body axes of the fore part of the fuselage denoted 
by xjyfZf. The component body axes can be obtained from XfyfZf through a sequence of rota- 
tions. For example, axes xfj/f-zf for the right half-wing can be obtained through a rotation 71 
about Xfy/Zf to an intermediate set of axes x' w y' w z' w and a rotation 72 about y' w to xf */f ~f , where 
in the case of the wing 72 is known as the dihedral angle. Table 2 gives the rotation angles for the 
individual components. 


Table 2 - Rotation Angles for Component Body Axes 


Component 

Rotations | 

Body Axes 

71 

72 

73 

xfyfzf 
Aft fuselage 

180° 

0 

0 

rpR^.R ~rR 

Right half-wing 

90° 

4° 

0 

Left half-wing 

-90° 

4° 

0 

Right horizontal stabilizer 

94.3CP 

9° 

0 

Left horizontal stabilizer 

-94.30P 

9° 

0 

V V V 
Ve z e 

Vertical Stabilizer 

0 

118.74? 

-90° 


Using analogies with the first of Eqs. (2), the various matrices of direction cosines are as follows: 


‘ -1 

0 

0 " 

II 

“o 3 

0 

0.9976 

-0.0698 “ 

0 

-1 

0 

-1 

0 

0 

0 

0 

1 

0 

0.0698 

0.9976 


ci 

ci 


0 -0.9976 -0.0698 

1 0 0 

0 -0.0698 0.9976 



-0.0750 0.9849 -0.1560 

-0.9972 -0.0741 0.0117 

0 0.1564 0.9877 


-0.0750 -0.9849 -0.1560 

0.9972 -0.0741 -0.0117 

0 -0.1564 0.9877 



-0.4808 0 -0.8768 

-0.8768 0 0.4808 

0 1 0 


(A2) 
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The inertia properties of the aircraft model are given in lumped form. To this end, the flexible 
components are divided into certain numbers of bays and the mass corresponding to each bay is 
lumped at the mass center of the bay. Moreover, the manner in which the mass is distributed over 
the associated cross-sectional area is represented by mass moments and mass products of inertia 
about axes with the origin at the mass center and parallel to the body axes of the respective com- 
ponent. Table 3 lists the coordinates of the mass centers, the mass values and the mass moments 
and mass products of inertia, in which the symmetry of the inertia matrices is implied. The masses 
have units lb • s 2 /in and the mass moments and mass products of inertia have units lb • s 2 • in. 
The cantilever beams lengths are L J = 295.86 in, Lf = 279.79 in, L % — L = 328.83 in, 
Lf = L^ = 127.46 in, L v e = 113.48 in. 


Table 3 - Inertia Properties 


Fore Part of the Fuselage 

No. 

'j' F 

x f 

yf 

z f 

m F 

T f 

J xxf 

T f 

J yyf 

T f 

J zzf 

7 f 

j *jlL 

jF 
J xz f 

jLi_ 

1 

7.41 

0 

5.20 

0.7237 

514.6542 

230.8500 

313.3728 

-0.1010 

-7.9229 

0 

2 

24.74 

0 

1.98 

0.2095 

206.4986 

115.1881 

97.2035 

-0.3495 

1.4603 

0 

3 

38.95 

0 

-1.54 

0.3014 

333.7084 

182.4735 

160.9599 

0.1424 

-1.9160 

0 

4 

48.49 

0 

3.25 

0.7558 

354.9372 

129.9051 

251.8536 

1.2350 

-0.6913 

0 

5 

55.08 

0 

8.60 

0.3555 

351.5220 

181.4974 

177.6291 

1.3775 

3.0682 

0 

6 

69.88 

0 

7.88 

0.5562 

559.9651 

274.9544 

306.9516 

-2.2111 

2.9413 

0 

7 

84.47 

0 

7.32 

0.3167 

298.6947 

182.3156 

124.5506 

-0.3573 

-2.2086 

0 

8 

101.30 

0 

5.29 

0.7053 

501.1153 

418.5278 

480.2490 

0.1657 

5.3648 

0 

9 

111.10 

0 

5.40 

0.1546 

137.1911 

84.7780 

54.2436 

-0.1605 

0.3754 

0 

10 

129.67 

0 

1.62 

1.0289 

934.8339 

556.8193 

535.0830 

-10.4733 

1.2842 

0 

11 

160.23 

0 

6.60 

1.2664 

919.9719 

477.2896 

595.6675 

3.1044 

17.3295 

0 

12 

189.38 

0 

-0.30 

1.0170 

673.8353 

340.9581 

437.3046 

2.3821 

-20.5815 

0 

13 

204.35 

0 

1.98 

1.0789 

623.3512 

354.2640 

386.3648 

-0.4686 

-1.1859 

0 

14 

223.81 

0 

1.16 

1.0497 

394.7331 

245.3210 

192.5843 

-0.8907 

9.7354 

0 

15 

237.71 

0 

6.55 

0.2496 

86.3807 

41.6213 

48.6587 

-0.5308 

0.0829 

0 

16 

251.30 

0 

6.57 

0.7234 

135.0187 

96.6805 

69.4034 

0.9425 

-9.6525 

0 

17 

266.50 

0 

14.02 

0.4702 

65.4238 

49.4018 

57.8918 

1.1289 

7.6562 

0 

18 

281.47 

0 

17.71 

0.0979 

14.8309 

8.2310 

10.2739 

-0.0104 

-0.0647 

0 

19 

292.08 

0 

17.66 

0.0640 

3.0164 

2.1283 

2.1801 

-0.0673 

0.0104 

0 

Aft Part of the Fuselage | 

No. 


yf 

z f 

A. 
m f 

J A 

Jxxf 

t a 

J yyf 

t a 

J zzf 

1 A 

J *vf 

j A 

J xzf 

J A 

-Jjy.r L- 

1 

1.67 

0 

0.81 

0.1434 

160.7476 

76.2673 

84.8272 

-0.0052 

0.1320 

0 

2 

10.66 

0 

17.72 

0.6147 

2114.0956 

253.7152 

1887.3572 

0.7172 

8.8913 

0 

3 

22.18 

0 

19.00 

1.4282 

3571.1313 

848.8131 

2879.2026 

7.9514 

-102.4285 

0 

4 

40.41 

0 

2.00 

0.1183 

133.0820 

72.5725 

60.7063 

0.0233 

0.0777 

0 

5 

46.82 

0 

8.04 

0.5808 

354.6912 

221.2985 

163.5387 

-0.1709 

-4.8263 

0 

6 

66.11 

0 

9.51 

0.2506 

259.3829 

186.0932 

84.9230 

0.4246 

-0.7664 

0 

7 

83.74 

0 

8.39 

0.7843 

900.7860 

470.6690 

478.1052 

8.3087 

-14.4348 

0 

8 

98.21 

0 

3.87 

0.5821 

317.6683 

222.4144 

138.4934 

-3.6275 

2.4908 

0 

9 

117.50 

0 

-8.33 

0.5870 

634.8913 

217.1299 

530.1143 

-0.3728 

11.5375 

0 

10 

142.62 

0 

-6.45 

0.6048 

317.3757 

146.8824 

212.2078 

0.5593 

-4.3913 

0 

11 

156.28 

0 

-10.77 

0.4945 

168.9475 

127.7664 

71.1511 

3.7414 

-6.0380 

0 

12 

173.22 

0 

-11.04 

0.1165 

62.1588 

48.9953 

21.8166 

0.4609 

1.1496 

0 

13 

189.26 

0 

-12.24 

0.2131 

49.4821 

37.8178 

19.5795 

0.5644 

1.3852 

0 

14 

203.41 

0 

-6.82 

0.1220 

36.0701 

28.0876 

15.0225 

-0.1605 

-3.7207 

0 

15 

218.31 

0 

-19.66 

0.0880 

40.5364 

37.9395 

11.4727 

-0.0129 

3.5368 

0 

16 

238.30 

0 

-10.17 

0.0585 

14.7040 

11.5064 

5.9785 

0.0078 

-0.8596 

0 

17 

253.16 

0 

-14.94 

0.0580 

6.6491 

6.2555 

3.7440 

0.1036 

-0.4220 

0 

18 

278.12 

0 

-17.35 

0.0098 

0.4143 

0.9580 

0.8415 

-0.0052 

-0.1087 

0 
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Table 3 - Continued 


Wing 

No. 


Vw 

z R 

___ ft 
TTX^ 

J R 

u xxw 

J R 

yyw 

J R 

u zzw 

J R 

xyw 

J R 

u xzw 

J R 

^ yzw 

1 

17.34 

-5.67 

-0.41 

0.7932 

592.7715 

135.2155 

675.5300 

-11.4074 

-20.6970 

23.8870 

2 

41.44 

-0.36 

0.64 

0.3380 

259.1733 

13.6761 

253.5778 

2.9790 

0.1402 

14.4486 

3 

55.80 

7.21 

0.72 

0.5481 

203.2508 

22.5351 

202.0929 

-2.4412 

-0.1781 

16.5185 

4 

68.45 

0.92 

-1.55 

0.2900 

145.9962 

10.7568 

144.9995 

1.1106 

0.4620 

4.8288 

5 

80.75 

4.29 

-0.63 

0.3517 

109.3290 

19.5652 

117.3550 

0.4698 

-1.1316 

4.4349 

6 

95.11 

2.34 

-0.89 

0.1960 

87.6687 

7.3973 

86.7071 

-3.5075 

-0.1369 

3.3911 

7 

109.21 

5.70 

-0.49 

0.1405 

68.3377 

5.1512 

67.9640 

0.3682 

-0.0707 

3.2896 

8 

123.58 

10.93 

-0.20 

0.2014 

98.0351 

8.6246 

100.8629 

-5.1787 

-0.2981 

4.7601 

9 

141.58 

6.17 

-0.36 

0.1441 

62.7718 

6.3306 

64.6489 

-0.3631 

-0.0224 

2.5948 

10 

159.69 

5.62 

-0.59 

0.1550 

58.9409 

7.3054 

62.3486 

-2.4040 

-0.0168 

1.5852 

11 

180.67 

2.46 

-0.60 

0.1481 

57.8061 

6.8510 

61.6321 

-1.7038 

0.0205 

1.7613 

12 

199.99 

5.62 

-0.49 

0.1504 

44.7968 

8.9780 

51.3853 

-4.0862 

-0.1240 

1.2184 

13 

223.12 

3.28 

-0.46 

0.1146 

26.2412 

4.9868 

29.6906 

0.2463 

0.0071 

0.5415 

14 

242.88 

1.87 

-0.49 

0.0825 

16.9014 

4.1168 

19.9998 

-0.2915 

-0.0464 

0.3456 

15 

260.53 

-0.02 

-0.71 

0.0700 

11.9896 

2.0960 

13.3314 

0.1592 

0.0681 

0.3213 

16 

278.05 

1.04 

-0.46 

0.0628 

9.2044 

1.9160 

10.6262 

0.1316 

0.0019 

0.0663 

17 

295.21 

-0.06 

-0.38 

0.0471 

6.9537 

1.6338 

8.2733 

0.2092 

-0.0066 

-0.0035 

18 

311.97 

-3.95 

-0.38 

0.0334 

4.6844 

0.6266 

5.1908 

0.1694 

-0.0063 

-0.0284 

19 

328.86 

0.04 

-0.40 

0.0287 

2.3874 

0.7263 

3.0305 

-0.0109 

-0.0045 

-0.0293 

| Horizontal Stabilizer | 

No. 

x e 

yf 


i? 
771 e 

Jxxe 

J R 

Jyye 

J R 

u zze 

Jxye 

J R 

u xze 

J R 

u yze 

1 

3.4154 

-2.0199 

4.8360 

0.1298 

152.7019 

24.5226 

24.9826 

9.4614 

-20.3966 

0.6193 

2 

14.6498 

4.1648 

-0.0644 

0.0618 

13.3528 

2.0288 

14.9172 

1.0676 

-0.0060 

0.0244 

3 

32.5615 

6.2772 

-0.0954 

0.0625 

12.1552 

1.9324 

13.7307 

0.2775 

0.0112 

0.0409 

4 

50.7350 

3.5463 

-0.1694 

0.0574 

8.3296 

2.7220 

10.7676 

0.5878 

0.0612 

-0.0078 

5 

66.0100 

3.8814 

-0.0389 

0.0273 

3.6287 

0.3475 

3.8652 

0.2049 

-0.0006 

0.0143 

6 

76.5595 

2.0048 

0.0296 

0.0173 

2.3491 

0.2014 

2.4673 

0.1894 

0.0006 

0.0055 

7 

88.2275 

2.7918 

-0.1213 

0.0348 

3.6527 

0.9998 

4.5598 

0.3189 

0.0164 

0.0305 

8 

105.1312 

2.4328 

-0.0898 

0.0259 

2.2711 

0.6068 

2.8183 

-0.0003 

-0.0048 

0.0187 

9 

121.3090 

3.3819 

-0.0468 

0.0199 

1.6192 

0.4246 

2.0069 

0.1574 

-0.0007 

0.0040 

10 

128.5213 

-2.0644 

-0.0276 

0.0135 

0.3576 

0.0114 

0.3586 

0.0197 

-0.0001 

-0.0010 

Vertical Stabilizer ] 

No. 

V 

vY 

zY 

V' 

m e 

J v 

°xxe 

J v 

Jyye 

J v 

u zze 

J v 

u xye 

J v 

J xze 

j v 

Jyze 

1 

-3.0130 

-4.0388 

0 

0.0867 

26.0215 

26.8758 

50.3184 

18.1691 

0.1131 

0.1417 

2 

17.3464 

-5.0753 

0 

0.0419 

8.2448 

2.9690 

10.2273 

3.9298 

0.0052 

0.0149 

3 

32.7117 

-3.0312 

0 

0.0347 

6.6885 

2.6766 

8.6842 

3.2180 

0.0097 

0.0124 

4 

48.0499 

-0.6415 

0 

0.0298 

7.0021 

2.2154 

8.7929 

3.2766 

0.0395 

0.0935 

5 

49.3577 

25.8074 

0 

0.0741 

5.0234 

42.2968 

46.9395 

-2.2970 

-0.0025 

-0.0045 

6 

90.2614 

-3.2365 

0 

0.0176 

1.7168 

1.1158 

2.6177 

1.0828 

-0.1246 

-0.1088 

7 

101.9491 

19.9499 

0 

0.0062 

0.2111 

0.1462 

0.3470 

0.0991 

0 

0 

8 

103.6312 

-3.8959 

0 

0.0150 

1.3610 

0.8632 

2.1464 

0.8286 

-0.0010 

0.0035 

9 

114.7248 

-10.6405 

0 

0.0065 

0.3567 

0.2077 

0.5619 

0.2227 

0 

0 

10 

116.1657 

2.4489 

0 

0.0523 

4.5419 

1.0689 

5.5228 

1.0377 

0 

0 

Engine Pylon (Attached to the Aft Part of the Fuselage) 

No. 

x f 

yf 

z f 

mf 

J A 

u xxf 

r A 

u yyf 

j A 

J zzf 

7 a 

J xyf 

T a 

J xzf 

T a 

J yzf 

1 

102.44 

-64.05 

-19.15 

3.0115 

347.7133 

1947.0343 

1872.5444 

-18.0506 

34.9386 

4.6139 

2 

108.62 

-37.00 

-13.96 

0.2647 

35.2545 

137.1807 

162.3658 

7.4957 

11.3756 

5.7131 


The stiffness properties consist of the flexural rigidity and torsional rigidity of the cross-sectional 
area at certain locations of the elastic components. Both have units lb • in 2 . Table 4 gives the loca- 
tions of the cross-sectional areas and the values of the corresponding rigidities. 
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Table 4 - Stiffness Properties 


| Fuselage | 

| Fore Part 

Aft Part 

No. 

X f 

EI fy „ EI u a 

Gjf 

No. 

x f 

EI fy ,, Elf s 

GJf 

1 

0.00 

1.10 x 10 11 7.40 x 10 1U 

4.07 x 10 10 

i 

1.29 

1.1 x 10“ 7.44 x lO 10 

4.1 x 10 1U 

2 

10.21 

1.10 x 10" 7.44 x 10 10 

4.10 x 10 10 

2 

5.04 

1.09 x 10" 7.40 x lO 10 

4.07 x lO 10 

3 

25.21 

1.10 x 10 11 7.47 x 10 10 

4.12 x 10 10 

3 

8.79 

1.04 x 10“ 7.37 x lO 10 

4.04 x lO 10 

4 

38.21 

1.10 x 10 u 7.47 x 10 10 

4.12 x 10 10 

4 

26.29 

9.44 x lO 10 6.87 x lO 10 

3.67 x lO 10 

5 

47.21 

1.10 x 10 11 7.47 x 10 10 

4.12 x 10 10 

5 

40.29 

7.98 x lO 10 6.17 x lO 10 

3.23 x lO 10 

6 

55.21 

1.10 x 10 11 7.47 x 10 10 

4.12 x 10 10 

6 

49.79 

6.72 x lO 10 5.32 x lO 10 

2.87 x lO 10 

7 

70.21 

1.10 x 10 11 7.45 x 10 10 

4.12 x lO 10 

7 

65.79 

6.72 x lO 10 4.15 x lO 10 

2.27 x lO 10 

8 

84.71 

1.10 x 10 11 7.36 x 10 10 

4.12 x 10 10 

8 

82.79 

6.45 x lO 10 3.21 x lO 10 

1.72 x lO 10 

9 

99.71 

1.10 x 10“ 7.20 x 10 10 

4.12 x 10 10 

9 

101.29 

5.15 x lO 10 2.43 x lO 10 

1.37 x lO 10 

10 

111.21 

1.10 x 10" 7.03 x 10 10 

4.12 x 10 10 

10 

120.29 

4.07 x lO 10 1.94 x lO 10 

1.12 x lO 10 

11 

129.21 

1.08 x 10" 6.80 x 10 l ° 

4.11 x 10 l ° 

11 

138.79 

2.30 x lO 10 1.37 x lO 10 

8.39 x 10 9 

12 

144.21 

1.02 x 10“ 6.54 x 10 l ° 

3.87 x 10 10 

12 

156.29 

1.09 x lO 10 9.68 x 10 9 

5.79 x 10 9 

13 

160.21 

9.06 x 10 10 6.14 x 10 10 

3.14 x 10 10 

13 

172.79 

1.63 x 10 l ° 6.62 x 10 9 

3.92 x 10 9 

14 

185.71 

6.5 x 10 10 5.25 x 10 10 

1.71 x lO 10 

14 

188.79 

2.10 x lO 10 4.39 x 10 9 

2.58 x 10 9 

15 

205.71 

4.37 x 10 10 4.23 x 10 l ° 

8.67 x 10 9 

15 

204.79 

1.38 x lO 10 2.80 x 10 9 

1.63 x 10 9 

16 

224.71 

2.68 x 10 10 3.97 x 10 10 

2.98 x 10 9 

16 

220.79 

7.48 x 10 9 1.70 x 10 9 

9.67 x 10® 

17 

238.11 

1.76 x 10 10 2.04 x 10 10 

2.19 x 10 9 

17 

237.29 

4.67 x 10 9 9.22 x 10® 

5.46 x 10® 

18 

251.11 

1.11 x 10 10 1.28 x 10 10 

1.27 x 10 9 

18 

256.29 

3.02 x 10 9 4.88 x 10® 

3.08 x 10® 

19 

267.06 

6.01 x 10 9 6.42 x 10 9 

6.43 x 10® 

19 

279.79 

1.30 x 10 9 2.40 x 10® 

1.17 x 10® 

20 

282.06 

3.13 x 10 9 2.87 x 10 9 

3.05 x 10® 





21 

295.86 

1.16 x 10 9 6.15 x 10 s 

7.72 x 10 7 





Wing 

| Horizontal Stabilizer 

No. 

A 

x w 

El* 

GJf 

No. 

x e 

Elf 

GJf 

1 

0.00 

1.09 x 10 l ° 

1.07 x 10 1(J 

1 

0.00 

3.92 x 10® 

2.43 x 10® 

2 

17.10 

9.70 x 10 9 

1.07 x lO 10 

2 

2.41 

3.78 x 10® 

2.36 x 10® 

3 

34.19 

8.65 x 10 9 

1.04 x lO 10 

3 

13.84 

3.16 x 10® 

2.09 x 10® 

4 

40.95 

8.10 x 10 9 

9.95 x 10 9 

4 

31.89 

2.32 x 10® 

1.73 x 10® 

5 

54.45 

6.95 x 10 9 

8.80 x 10 9 

5 

49.94 

1.66 x 10® 

1.35 x 10® 

6 

67.85 

5.70 x 10 9 

7.60 x 10 9 

6 

64.33 

1.23 x 10® 

1.03 x 10® 

7 

81.11 

4.84 x 10 9 

6.40 x 10 9 

7 

75.11 

9.70 x 10 7 

8.05 x 10 7 

8 

94.45 

4.28 x 10 9 

5.30 x 10 9 

8 

88.35 

7.10 x 10 7 

5.85 x 10 7 

9 

107.95 

3.49 x 10 9 

4.50 x 10 9 

9 

103.98 

4.82 x 10 7 

4.14 x 10 7 

10 

122.95 

2.81 x 10® 

3.92 x 10 9 

10 

119.63 

3.27 x 10 7 

2.98 x 10 7 

11 

140.45 

2.18 x 10 9 

3.19 x 10 9 

11 

127.46 

2.74 x 10 7 

2.19 x 10 7 

12 

159.45 

1.68 x 10 9 

2.43 x 10 9 

Vertical Stabilizer 

13 

179.45 

1.35 x 10 9 

1.86 x 10 9 

No. 

xX 

Elf 

GJf 

14 

200.70 

1.02 x 10 9 

1.37 x 10 9 

1 

0.00 

1.56 x 10 9 

7.30 x 10® 

15 

221.95 

7.35 x 10® 

8.80 x 10® 

2 

7.40 

1.33 x 10 9 

6.70 x 10® 

16 

241.70 

5.65 x 10® 

5.55 x 10® 

3 

21.15 

9.66 x 10® 

5.57 x 10® 

17 

259.95 

4.28 x 10® 

3.56 x 10® 

4 

34.90 

6.86 x 10® 

4.45 x 10® 

18 

276.70 

3.54 x 10® 

2.53 x 10® 

5 

49.73 

4.61 x 10® 

3.55 x 10® 

19 

293.95 

2.82 x 10® 

1.71 x 10® 

6 

63.98 

3.06 x 10® 

2.41 x 10® 

20 

311.10 

2.02 x 10® 

8.55 x 10 7 

7 

77.63 

2.01 x 10® 

1.38 x 10® 

21 

327.01 

1.27 x 10® 

6.50 x 10® 

8 

91.31 

1.26 x 10® 

8.69 x 10 7 





9 

105.82 

6.88 x 10 7 

3.56 x 10 7 





10 

113.48 

4.25 x 10 7 

8.51 x 10® 


The aerodynamic forces involve the slope C La , of the lift curve, the drag coefficient C Di0 , 
the slope C S Qi of the lateral force curve and the control effectiveness coefficients C L s a , C s s e and 
C s sr- To determine the aerodynamic forces, the aircraft components are divided into a given num- 
ber of sections and the coefficients are given for each of the sections. A typical section has a 
trapezoidal shape defined by two points x ia y ia z ia and XibVibZib and by respective chords c, a and 
Cib, where the chords are parallel to the longitudinal axis xj of the fuselage (Fig. A). For all 
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sections, Cp/o = Cdwo = 0.016 and kf = k w = 0.04. The lines of the aerodynamic centers are 
also shown in Fig. A; they are located at one quarter of the chord of the wing and horizontal and 
vertical stabilizers. The line of aerodynamic centers for the fuselage are located at one half of the 
chord of each aerodynamic section, as shown in Fig. A. The aerodynamic forces for the fore and 
aft parts of the fuselage are acting at the line of the aerodynamic centers. The control forces, on 
the other hand, are acting at 0.55 of the chord. The aerodynamic coefficients corresponding to 
the sections just mentioned are listed in Table 5. The slope of the lift coefficients and the control 
effectiveness for the aileron and the elevator listed in Table 5 are given for a single Mach number 
of 0.75. They must be corrected for different Mach numbers by the compressibility factor given in 
Fig. A. Note that there are only seven sections for the wing with control effectiveness coefficients; 
they correspond to the aileron. 


Table 5 




Fore Part of the Fuselage - 

Horizontal Lifting Surface 



Xfa 

F 
V fa 

*fa 

- c /o - 

Xfb 

Vfb 

7 F 

_Jjb_ 

c n b 

r F 

^ Laf 

111.21 

-24.00 

0.17 

40.00 

111.21 

24.00 

0.17 

40.00 

0.065 

151.21 

-36.25 

0.17 

40.00 

151.21 

36.25 

0.17 

40.00 

0.065 

191.21 

-36.25 

0.17 

40.00 

191.21 

36.25 

0.17 

40.00 

0.20 

231.21 

-36.25 

0.17 

40.00 

231.21 

36.25 

0.17 

40.00 

0.50 

277.21 

-36.25 

0.17 

46.00 

277.21 

36.25 

0.17 

46.00 

0.50 



Fore Part of the Fuselage 

- Vertical Lifting Surface 


j 

X fQ 

Vfa 

z fa 


X fb 

Vfb 


F 

C fb 

c Ui 

13.21 

0 

41.17 

56.00 

13.21 

0 

41.17 

56.00 

0.142 

69.21 

0 

39.17 

56.00 

69.21 

0 

39.17 

56.00 

0.142 

126.23 

0 

37.17 

56.80 

126.23 

0 

37.17 

56.80 

0.142 

182.21 

0 

35.17 

55.98 

182.21 

0 

35.17 

55.98 

0.242 

248.21 

0 

33.17 

57.00 

248.21 

0 

33.17 

57.00 

0.305 

Wing | 

x wa 

Vwa 

^wa 

c wa 

x wb 

Vwb 

^R 

wb 

Cwb 

ClOlUj 

CLSa 

0 

-62.49 

0 

129.49 

34.00 

-62.49 

0 

129.49 

4.9675 


34.00 

-62.49 

0 

129.49 

58.36 

-52.40 

0 

115.51 

6.1879 


58.36 

-52.40 

0 

115.51 

80.71 

-43.15 

0 

102.68 

6.5604 


80.71 

-43.15 

0 

102.68 

101.07 

-34.72 

0 

90.99 

6.7609 


101.07 

-34.72 

0 

90.99 

120.35 

-33.06 

0 

86.25 

6.9328 


120.35 

-33.06 

0 

86.25 

138.62 

-31.49 

0 

81.75 

7.1620 


138.62 

-31.49 

0 

81.75 

155.90 

-30.00 

0 

77.50 

7.4485 


155.90 

-30.00 

0 

77.50 

172.17 

-28.60 

0 

73.50 

7.7349 


172.17 

-28.60 

0 

73.50 

187.43 

-27.29 

0 

69.75 

7.9641 


187.43 

-27.29 

0 

69.75 

203.59 

-25.89 

0 

65.77 

8.1647 

2.4494 

203.59 

-25.89 

0 

65.77 

219.20 

-24.55 

0 

61.93 

8.2792 

2.4838 

219.20 

-24.55 

0 

61.93 

234.30 

-23.25 

0 

58.21 

8.3365 

2.5010 

234.30 

-23.25 

0 

58.21 

248.91 

-21.99 

0 

54.62 

8.3365 

2.5010 

248.91 

-21.99 

0 

54.62 

263.01 

-20.78 

0 

51.15 

8.2506 

2.4752 

263.01 

-20.78 

0 

51.15 

276.62 

-19.61 

0 

47.80 

8.0214 

2.4064 

276.62 

-19.61 

0 

47.80 

289.72 

-18.48 

0 

44.58 

7.7063 

2.3119 

289.72 

-18.48 

0 

44.58 

301.25 

-17.49 

0 

41.74 

7.2766 


301.25 

-17.49 

0 

41.74 

311.27 

-16.62 

0 

39.28 

6.6463 


311.27 

-16.62 

0 

39.28 

320.29 

-15.85 

0 

37.06 

5.7869 


320.29 

-15.85 

0 

37.06 

327.66 

-15.21 

0 

35.25 

4.5837 


327.66 

-15.21 

0 

35.25 

335.02 

-14.58 

0 1 

33.43 

2.5783 
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Table 5 (Continued) 
Horizontal Stabilizer 


^ea 

Vea 

^ ea 

u ea 

E 

Xeb 

Veb 

z eb 

C e b 

C*Lae 

CLSe 

0 

-25.40 


65.31 



0 

61.62 

1.5783 

0.7891 

10.98 

-24.06 


61.92 

23.37 

-22.77 

0 

58.67 

1.6844 

0.8422 

23.37 

-22.77 


58.67 

35.24 

-21.54 

0 

55.55 

1.8104 

0.9052 

35.24 

-21.54 


55.55 



0 

52.57 

1.9828 

0.9414 

46.60 

-20.36 


52.57 

57.44 

-19.23 

0 

49.72 

2.2282 

1.1141 

57.44 

-19.23 


49.72 

67.88 

-18.14 

0 

46.98 

2.5398 

1.2699 

67.88 

-18.14 


46.98 

77.91 


0 

44.34 

2.8913 

1.4457 

77.91 

-17.10 


44.34 

87.53 


0 

41.82 

3.1632 

1.5816 

87.53 

-16.10 


41.82 

96.74 

-15.15 

0 


3.3025 

1.6512 

96.74 

-15.15 


■ 


-14.23 

0 


3.2958 

1.6479 

105.54 

-14.23 


■ 

113.93 

-13.36 

0 

34.88 



113.93 

-13.36 


34.88 

121.91 

-12.53 

0 

32.79 

2.7454 

1.3727 

121.91 

-12.53 


32.79 


-11.79 

0 


1.9828 


Vertical Stabilizer ] 

x v 

x eo 

Via 


c ea 

x eb 

v 

Veb 

•4 

V 

^eb 

r v 

' s&e 

n v 

s6r 

-31.15 

-39.58 

0 

102.00 

-20.81 

-37.96 

0 

97.85 

0.9482 

0.4267 

-20.81 

-37.96 

0 

97.85 

-4.81 

-35.44 

0 

91.43 

1.4740 

0.6633 

-4.81 

-35.44 

0 

91.43 

3.05 

-34.20 

0 

88.28 

1.7290 

0.7781 

3.05 

-34.20 

0 

88.28 

8.92 

-33.28 

0 

85.92 

2.6550 

1.1948 

8.92 

-33.28 

0 

85.92 

25.14 

-30.73 

0 

79.42 

2.8540 

1.2843 

25.14 

-30.73 

0 

79.42 

41.37 

-28.18 

0 

72.91 

3.1030 

1.3964 

41.37 

-28.18 

0 

72.91 

56.34 

-25.82 

0 

66.90 

3.3020 

1.4859 

56.34 

-25.82 

0 

66.90 

71.32 

-23.47 

0 

60.89 

3.4280 

1.5426 

71.32 

-23.47 

0 

60.89 

86.30 

-21.11 

0 

54.88 

3.4350 

1.5458 

86.30 

-21.11 

0 

54.88 

103.37 

-18.43 

0 

48.03 

3.2020 

1.4409 

103.37 

-18.43 

0 

48.03 

113.44 

-16.84 

0 

43.99 

2.3430 

1.0544 

Engine Pylon 

a 

Xfa 

yfa 

z fa 

C /a 

x fb 

yfb 

A 

z fb 

c % 

°Lq/ 

50.79 

0.00 

-5.13 

126.00 

50.79 

20.00 

-11.83 

126.00 

1.2376 

50.79 

20.00 

-11.83 

126.00 

50.79 

33.00 

-16.16 

117.33 

1.2892 

50.79 

33.00 

-16.16 

117.33 

50.79 

47.00 

-20.83 

108.00 

1.2900 

42.79 

47.00 

-20.83 

116.00 

42.79 

56.00 

-23.84 

116.00 

1.1250 

42.79 

56.00 

-23.84 

116.00 

42.79 

65.00 

-26.85 

116.00 

1.0500 

42.79 

65.00 

-26.85 

116.00 

42.79 

74.00 

-29.86 

116.00 

0.8500 

42.79 

74.00 

-29.86 

116.00 

42.79 

82.00 

-32.53 

116.00 

0.5500 


69 








REPORT DOCUMENTATION PAGE 


Form Approved 



Public reporting burden lor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data 
sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations 
and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704- 
0188), Washington, DC 20503. 

1 . AGENCY USE ONLY ( Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

July 2003 Contractor Report 

4. TITLE AND SUBTITLE 

Integrated Approach to the Dynamics and Control of Maneuvering Flexible 
Aircraft 

5. FUNDING NUMBERS 

NCC 1-346 
706-31-41-02 

6. AOTHOR(S) 

Leonard Meirovitch and Ilhan Tuzcu 

s 7. PERM)RMIN6 6 r6aNIZATION NAME(S) AND ADDRESS(E§) “ “ 

Virginia Polytechnic Institute and State University 
Engineering Science and Mechanics Department, MC 0219 
Blacksburg, VA 24061 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 2368 1-2199 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

N ASA/CR-2003-2 1 1748 


11. SUPPLEMENTARY NOTES 

An electronic version can be found at http://techreports.larc.nasa.gov/ltrs/ or http://ntrs.nasa.gov 
Langley Technical Monitor: Martin R. Waszak 


12a. DISTRIBUTION/AVAILABILITY STATEMENT I 12b. DISTRIBUTION CODE 


Unclassified-Unlimited 

Subject Category 08 Distribution: Standard 

Availability: NASA CASI (301) 621-0390 


13. ABSTRACT (Maximum 200 words) 

This work uses a fundamental approach to the problem of simulating the flight of flexible aircraft. To this end, it 
integrates into a single formulation the pertinent disciplines, namely, analytical dynamics, structural dynamics, 
aerodynamics and controls. It considers both the rigid body motions of the aircraft, three translations (forward 
motion, sideslip and plunge) and three rotations (roll, pitch and yaw), and the elastic deformations of every point 
of the aircraft, as well as the aerodynamic, propulsion, gravity and control forces. The equations of motion are 
expressed in a form ideally suited for computer processing. A perturbation approach yields a flight dynamics 
problem for the motions of a quasi-rigid aircraft and an "extended aeroservoelasticity" problem for the elastic 
deformations and perturbations in the rigid body motions, with the solution of the first problem entering as an 
input into the second problem. The control forces for the flight dynamics problem are obtained by an "inverse" 
process and the feedback controls for the extended aeroservoelasticity problem are determined by the LQG 
theory. A numerical example presents time simulations of rigid body perturbations and elastic deformations 
about 1) a steady level flight and 2) a level steady turn maneuver. 


14. SUBJECT TERMS 

Flexible Aircraft, Integrated Approach, Flight Dynamics, 
Extended Aeroservoelasticity, Time Simulations 


15. SECUR IT Y CLASSI F ICA T ION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION " 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

76 


16. PRICE CODE 


59. LIMIT A 1 I ON 

OF ABSTRACT 

UL 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


standard Form 298 (Rev. 2-89 

Prescribed by ANSI Std. Z-39-18 
298-102 


NSN 7540-01-280-5500 















